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Abstract 



In this Thesis we develop a novel direction in the theory of nano-objects, i.e., struc- 
tures of nanometer size in a tunnel contact with macroscopic electron reservoirs (metallic 
leads). This theory arose and was developed rapidly during the two recent decades as a 
response to challenging achievements of modern nanotechnology and experimental tech- 
niques. Evolution of this technology enabled the fabrication of various low-dimensional 
systems from semiconductor heterostructures to quantum wires and constrictions, quan- 
tum dots, molecular bridges and artificial structures constructed from large molecules. 
This impressive experimental progress initiated the development of a new direction in 
quantum physics, namely, the physics of artificial nano-objects. 

We focus in this work on a theoretical investigation of the Kondo physics in quantum 
dots and molecules with strong correlations. An exciting series of recent experiments on 
mesoscopic and nanoscale systems has enabled a thorough and controlled study of basic 
physical problems dealing with a local moment interacting with a Fermi sea of conduction 
electrons. Scanning tunnel microscopy and quantum-dot devices have provided new tools 
for studying the Kondo effect in many new perspectives and with unprecedented control. 
As experimental and theoretical investigations of tunneling phenomena continue, it turns 
out that the physics of tunneling spectroscopy of large molecules and complex quantum 
dots have much in common. 

In particular we elucidate the Kondo effect predicted in tunneling through triple quan- 
tum dots and sandwich-type molecules adsorbed on metallic substrate, which are referred 
to as trimers. The unusual dynamical symmetry of nano-objects is one of the most intrigu- 
ing problems, which arise in the theory of these systems. We demonstrate that trimers 
possess dynamical symmetries whose realization in Kondo tunneling is experimentally tan- 
gible. Such experimental tuning of dynamical symmetries is not possible in conventional 
Kondo scattering. We develop the general approach to the problem of dynamical symme- 
tries in Kondo tunneling through nano-objects and illustrate it by numerous examples of 
trimers in various configurations, in parallel, in series and in ring geometries. 

In the first part of this Thesis, the evenly occupied trimer in a parallel geometry is 
studied. We show that Kondo tunneling through the trimer is controlled by a family of 
SO{n) dynamical symmetries. The most striking feature of this result is that the value of 
the group index n = 3,4,5,7 can be changed experimentally by tuning the gate voltage 
applied to the trimer. Following the construction of the corresponding o„ algebras, the 
scaling equations are derived and the Kondo temperatures are calculated. 



In the second part, the Kondo physics of trimer with both even and odd electron 
occupation in a series geometry is discussed. We derive and solve the scaling equations 
for the evenly occupied trimer in the cases of the P x SO{A) x 5'0(4), -S'0(5), SO{l) and 
P X S'0(3) X S0{?)) dynamical symmetries. The dynamical-symmetry phase diagram is 
displayed and the experimental consequences are drawn. The map of Kondo temperature 
as a function of gate voltages is constructed. In addition, the influence of magnetic fleld 
on the dynamical symmetry and its role in the Kondo tunneling through the trimer are 
studied. It is shown that the anisotropic Kondo effect can be induced in the trimer by an 
external magnetic field. The corresponding symmetry group is SU (3). In the case of odd 
electron occupation, the effective spin Hamiltonian of the trimer manifests a two-channel 
Kondo problem albeit only in the weak coupling regime (due to unavoidable anisotropy) . 

In the third part, the point symmetry C^y of an artificial trimer in a ring geometry 
and its interplay with the spin rotation symmetry SU{2) are studied. This nano-object 
is a quantum dot analog of the Coqblin-Schrieffer model in which the Kondo physics is 
governed by a subtle interplay between spin and orbital degrees of freedom. The orbital 
degeneracy is tuned by a magnetic fleld, which affects the electron phases thereby leading 
to a peculiar Aharonov-Bohm effect. 

The following novel results were obtained in the course of this research: 

• It is found that evenly occupied trimer manifests a new type of Kondo effect that 
was not observed in conventional spin 1/2 quantum dots. The dynamical SO{n) 
symmetries of Kondo tunneling through evenly occupied trimer both in parallel and 
series geometry are unravelled. These symmetries can be experimentally realized 
and the specific value of n = 3, 4, 5, 7 can be controlled by gate voltage and/or tun- 
neling strength. The Kondo temperature explicitly depends on the index n and this 
dependence may be traced experimentally by means of measuring the variation of 
tunnel conductance as a function of gate voltages. The hidden dynamical symmetry 
manifests itself, firstly in the very existence of the Kondo effect in trimer with even 
occupation, secondly in non-universal behavior of the Kondo temperature Tk- In a 
singlet spin state the anisotropic Kondo effect can be induced in the trimer by an 
external magnetic field. 

• It is shown that the effective spin Hamiltonian of a trimer with odd electron occupa- 
tion weakly connected in series with left {I) and right (r) metal leads is composed of 

two-channel exchange and co-tunneling terms. Renormalization group equations for 
the corresponding three exchange constants J;, J,- and Jir are solved in a weak cou- 
pling limit (single loop approximation). Since J;^ is relevant, the system is mapped 
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on an anisotropic two-channel Kondo problem. The structure of the conductance 
as function of temperature and gate voltage implies that in the weak and interme- 
diate coupling regimes, two-channel Kondo physics persists at temperatures as low 
as several Tp^. Analysis of the Kondo effect in cases of higher spin degeneracy of 
the trimer ground state is carried out in relation with dynamical symmetries. The 
Kondo physics remains that of a fully screened impurity, and the corresponding 
Kondo temperature is calculated. 

• It is demonstrated that spin and orbital degrees of freedom interlace in ring shaped 
artificial trimer thereby establishing the analogy with the Coqblin-Schrieffer model 

in magnetically doped metals. The orbital degrees of freedom are tunable by an 
external magnetic field, and this implies a peculiar Aharonov-Bohm effect, since the 
electron phase is also affected. The conductance is calculated both in three- and 
two-terminal geometries. It is shown that it can be sharply enhanced or completely 
blocked at definite values of magnetic fiux through the triangular loop. 
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Chapter 1 
Introduction 



Background. Recently, studies of the physical properties of artificially fabricated nano- 
objects turn out to be a rapidly developing branch of fundamental and applied physics. 
Progress in these fields is stimulated both by the achievements of nanotechnology and 
by the ambitious projects of information processing, data storage, molecular electronics 
and spintronics. The corresponding technological evolution enabled the fabrication of 
various low-dimensional systems from semiconductor heterostructures to quantum wires 
and constrictions, quantum dots (QD), molecular bridges and artificial structures with 
large molecules built in electric circuits [1, 2J. This impressive experimental progress led to 
the development of nanophysics, a new aspect and research direction in quantum physics 
[3j. Artificial nano-objects possess the familiar features of quantum mechanical systems, 
but sometimes one may create in artificially fabricated systems such conditions, which 
are hardly observable "in natura". For example, one-dimensional to two-dimensional 
{ID 2D) crossover may be realized in quantum networks jH El E] and constrictions 
[7, 8j. The Kondo effect may be observed in non-equilibrium conditions [3 CHI at 
high magnetic fields [121 CHI El 113 CEl , and at finite frequencies [m UHlliniilll El 122] • 
Moreover, a quantum dot in the Kondo regime can be integrated into a circuit exhibiting 
the Aharonov-Bohm effect ^^1^1^ . 

According to the theory of Kondo effect in QD [221 12Z| , the spin degrees of the QD are 
involved in Kondo resonance, and Kondo effect takes place only if the dot has a nonzero 
spin in the ground state. Resonance Kondo tunneling was experimentally observed in 
QD with odd electron occupation number and spin one-half ground state fl ^ 1 ^ IHI H m]. 
But real nano-objects cannot be simply represented by a spin 1/2 moment, because low- 
lying spin excitations, which are always present in few-electron systems, should be taken 
into account. Evenly occupied dot usually has a singlet ground state. However, it was 
predicted theoretically ^12, and confirmed experimentally that Kondo tunneling can 
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be induced by an external magnetic field in planar QD with even occupation. This 
unconventional magnetic field induced Kondo effect arises because the spectrum of the 
dot possesses a low-lying triplet excitation above the singlet ground state. The Zeeman 
splitting energy of a triplet in an external magnetic field may exactly compensate the 
energy spacing between the two adjacent levels, and the lowest spin excitation possesses an 
effective spin 1/2, thus inducing a Kondo-like zero-bias anomaly (ZBA) in the differential 
conductance. A similar scenario may be realized in vertical QDs ITK| ITIH 122] where 
the Larmor (instead of the Zeeman) effect comes into play. The influence of an external 
magnetic fleld on the orbital part of the wave functions of electrons in vertical QDs is, 
in general, more pronounced than the Zeeman effect. Hence, singlet-triplet level crossing 
are induced by this magnetic fleld, causing the emergence of Kondo scattering 

Another device which manifests the Kondo effect in QDs with even electron occupa- 
tion is a double quantum dot. Double quantum dots (DQDs), namely, quantum dots with 
two potential-wells, oriented parallel to the lead surface were fabricated several years ago 
[SHI El]- The two wells in a DQD may be identical or have different size; the DQD may be 
integrated within an electric circuit either in series or in parallel; different gate voltages 
may be applied to each well. Moreover, one of the two wells may be disconnected from 
the leads (side geometry) EEl EII ■ If the tunneling between the right and left wells 
of the DQD is taken into account, DQD can be treated as an artificial molecule where 
the interdot tunneling results in the formation of complicated manifold of bonding and 
antibonding states which modifies its degrees of freedom |SH1 EHl SHI E] • The systematic 
treatment of the physics of DQD coupled to metallic leads [SHI EE] is based on the mecha- 
nism according to which the transition from a singlet state in a weak coupling regime to a 
triplet state in a strong coupling regime is an intrinsic property of nano-objects with even 
occupation. It is manifested in tunneling through real and artificial molecules in which 
the electrons are spatially separated into two groups with different degree of localiza- 
tion. Electrons in the first group are responsible for strong correlation effects (Coulomb 
blockade), whereas those in the second group are coupled to a metallic reservoir. The 
necessary precondition under which the singlet = ground state changes into a par- 
tially screened triplet S = 1 Kondo state due to hybridization with metallic leads is the 
existence of charge-transfer exciton in the DQD. Unlike quantum dot with odd occupa- 
tion whose Hamiltonian is mapped on the Kondo-type S(i-exchange Hamiltonian with a 
localized spin S = 1/2 obeying SU{2) symmetry, DQD in contact with metallic leads can 
be treated as a quantum spin rotator with S = 1. The effective Hamiltonian of DQD 
possesses the dynamical symmetry 5*0(4) of a spin rotator [HH] . 

The Kondo physics seems to be richer in systems involving tunneling through artificial 
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molecules containing more than two wells. Meanwhile, the Kondo effect was observed 
already in complex molecules containing cages with magnetic ions ^\ H3] . Analogy 
between the Kondo effect in real and artificial molecules in tunnel contact with metallic 
leads was noted some time ago in Refs. [33 EHl EHI- There exists a great variety of 
molecules containing magnetic rare earth (RE) ions secluded in carbon and nitrogen based 
cages. EndofuUerenes REC82 are the most common among them O ^1 |ini EZl ■ In these 
molecules magnetic ions are inserted in a nearly spherical carbon cage. Lanthanocenes 
Ln(C8H8)2 are sandwich-type molecules formed by two rings of CH radicals and magnetic 
ion Ln=Ce, Yb in between ^Hl UHl ED]- In these molecules the electrons in a strongly 
correlated / shell of Ln are coupled with loosely bound tt electrons in the cage. The ground 
state of cerocene molecule is a spin singlet combination ^Aig^fir^) of an / electron and 
TT orbitals, and the energy of the first excited triplet state ^E2g is rather small (~ 0.5 eV). 
In the ytterbocene (hole counterpart of cerocene) the ground state with one / hole is 
a triplet, and the gap for a singlet excitation is tiny, ~ 0.1 eV. In all these systems 
there is no direct overlap between the strongly correlated / electrons and the metallic 
reservoir. However, these electrons can influence the tunnel properties of the molecule 
via covalent bonding with outer vr electrons which are coupled to the metallic reservoir. 
Other examples of molecules secluding magnetic ions may be found, e.g., in HH]. 

Objectives. In this thesis we develop a theory of Kondo tunneling through triple 
quantum dots and sandwich-type molecules adsorbed on metallic substrate, which are 
referred to as trimers. One of the most intriguing phenomenon which arise in the theory 
of these systems is the emergence of an unusual dynamical symmetry of nano- objects. 
Our main purpose is to demonstrate that trimers possess dynamical symmetries whose 
realization in Kondo tunneling is experimentally tangible. Such experimental tuning of dy- 
namical symmetries is not possible in conventional Kondo scattering. In many cases even 
the very existence of Kondo tunneling crucially depends on the dynamical spin symmetry 
of trimer. We develop the general approach to the problem of dynamical symmetries in 
Kondo tunneling through nano-objects and illustrate it by numerous examples of trimer in 
various configurations, in parallel, in series and in ring geometries. We show that Kondo 
tunneling reveals hidden SO{n) dynamical symmetries of evenly occupied trimers both in 
parallel and series geometry. The possible values = 3, 4, 5, 7 can be controlled by gate 
voltages, indicating that abstract concepts such as dynamical symmetry groups are exper- 
imentally realizable. We construct the corresponding o„ algebras, derive and solve scaling 
equations and calculate the Kondo temperatures. We elucidate the role of discrete and 
continuous symmetries exposed by the Kondo effect in artificial trimer in ring geometry, 
i.e., triangular triple quantum dot (TTQD). In comparison with a linear configuration of 
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three quantum dots, a TTQD possesses additional degrees of freedom, namely, discrete 
rotations. We show that the Kondo physics of TTQD is determined by a subtle interplay 
between continuous spin rotation symmetry SU{2), and discrete point symmetry C^y. 
Moreover, such ring shaped nano-object can serve both as a Kondo-scatterer and as a 
peculiar Aharonov-Bohm (AB) interferometer, since the magnetic flux affects not only 
the electron phase, but also the nature of the ground and excited states of the trimer. 
The main lesson to be learned is that Kondo physics in trimer suggests a novel and in 
some sense rather appealing aspect of low-dimensional physics of interacting electrons. It 
substantiates, in a systematic way, that dynamical symmetry groups play an important 
role in mesoscopic physics. In particular, we encounter here some "famous" groups which 
appear in other branches of physics. Thus, the celebrated group SU{3) enters also here 
when a trimer is subject to an external magnetic field. And the group SO{5) which plays 
a role in the theory of superconductivity is found here when a certain tuning of the gate 
voltages in trimer is exercised. 

Structure. The structure of the Thesis is as follows. In the second Chapter, the basic 
physics of Kondo effect in quantum dots is brieffy reviewed. First, the Kondo tunneling 
through single quantum dot with odd electron occupation is described. Next, quantum 
dot with even number of electrons, subject to an external magnetic field is considered. 
It is shown that the system exhibits Kondo effect in a finite magnetic field, when the 
Zeeman energy is equal to the single-particle level spacing in the dot. Then, the Kondo 
physics of evenly occupied double quantum dot (DQD) is presented. Special attention 
is given to the symmetry properties of the DQD. It is shown that the DQD possesses 
the SO (4) dynamical symmetry of a spin rotator. Finally, we explain the concept of 
dynamical symmetry and its realization in complex quantum dots. 

In the third Chapter, the special case of trimer with even electron occupation in the 
parallel geometry is studied. We discuss the energy spectrum of the isolated trimer, derive 
renormalization group equations and demonstrate possible singlet-triplet level crossing 
due to tunneling. We show that the trimer manifests SO{n) dynamical symmetry in the 
Kondo tunneling regime. We expose the effective spin Hamiltonian of the trimer and 
construct the corresponding o„ algebras for the P x 5'0(4) x SO{A), SO{5) and SO{7) 
dynamical symmetries. We derive scaling equations and calculate the Kondo temperatures 
for the cases of P x 5'0(4) x 5'0(4) and SO{5) symmetries. 

In the fourth Chapter we discuss the physics of trimer in a series geometry and point 
out similarities and differences between Kondo physics in the parallel and series geome- 
tries. The scaling equations are derived and the Kondo temperatures are calculated for 
the evenly occupied trimer in the cases of the P x SO (4) x 5'0(4), SO (5), SO (7) and 



6 



P X SO{?)) X SO{?>) dynamical symmetries. The dynamical symmetries of the trimer are 
summarized by a phase diagram which can be scanned experimentally by appropriate vari- 
ations of gate voltages. We discuss a novel phenomena, namely, a Kondo effect without 
a localized spin. The anisotropic exchange interaction occurs between the metal electron 
spin and the trimer Runge-Lenz operator alone in an external magnetic field. The sym- 
metry group for such magnetic field induced anisotropic Kondo tunneling is S'f/(3). We 
show that in the case of odd occupation, the effective Hamiltonian of the trimer manifests 
generic futures of a two-channel Kondo problem at least in the weak coupling regime. 

In the fifth Chapter we concentrate on the point symmetry C^y of artificial trimer in 
a ring geometry and its interplay with the spin rotation symmetry SU{2) in the context 
of Kondo tunneling through triangular artificial molecule. The underlying Kondo physics 
is determined by the product of a discrete rotation symmetry group in real space and 
a continuous rotation symmetry in spin space. These symmetries are reflected in the 
resulting exchange Hamiltonian which naturally involves spin and orbital degrees of free- 
dom thereby establishing the analogy between the Coqblin-Schrieffer model in real metals 
and the physics of transport in complex quantum dots. The ensuing poor-man scaling 
equations are solved and the Kondo temperature is calculated. We show that the trimer 
also reveals a peculiar Aharonov-Bohm effect where the magnetic field affects not only 
the electron phase but also controls the underlying dynamical symmetry group. 

The derivation of the pertinent effective spin Hamiltonians and the establishment of 
group properties (in particular identification of the group generators and checking the cor- 
responding commutation relations) sometimes require lengthy mathematical expressions. 
These are collected in the appendices. 

This work was partially presented by posters and lectures in scientific conferences and 
schools (see List of Presentations). The first part of the results was published in Refs. 1, 
2 (see List of Publications). The second part was published in Refs. 3-5. The third part 
was published in Refs. 6, 7. 
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Chapter 2 

Kondo Physics in Artificial 
Nano-objects 



In this Chapter we review the basic physics of Kondo effect in single and double quantum 
dots. Quantum dot with odd electron occupation in tunnel contact with metallic leads is 
considered in Section \2.1\ The Hamiltonian of the system is written down within the 
framework of the Anderson model. The Haldane renormalization group procedure is de- 
scribed and the effective spin Hamiltonian is obtained by means of the Schrieffer- Wolff 
transformation. The expressions for the Kondo temperature and zero-bias conductance of 
the dot are then derived. In Section \2.2 the possibility of Kondo effect in evenly occupied 
single quantum dots subject to an external magnetic field is discussed. The Kondo effect 
in evenly occupied double quantum dot (DQD) is studied in Section \2. ,91 It is shown that 
DQD possesses S0{4) dynamical symmetry of a spin rotator. Finally in Section "KJ 



we 



introduce the concept of dynamical symmetry and its emergence in complex quantum dots. 



2.1 Kondo Effect in Single Quantum Dot (QD) 

The conventional Kondo effect appears in scattering of conduction electrons by localized 
magnetic impurity 51J. The latter is represented by its spin S, which possesses the 
SU{2) symmetry of rotationally invariant moment. Spin scattering of conduction electrons 
dynamically screens this moment, and the system transforms into a local Fermi liquid 
with separate branches of charge and spin excitations ES] • According to Refs. [23 
n7\ . the problem of tunneling through a nano-object with odd electron occupation and 
strong Coulomb blockade suppressing charge fluctuations can be mapped on the Kondo 
scattering problem. The Kondo effect emerges in a quantum dot occupied by an odd 
number of electrons at temperatures below the mean level spacing in the dot. Under such 



10 



conditions, the highest occupied level Eq filled by a single electron produces the Kondo 
effect. The other levels, occupied by pairs of electrons with opposite spins, don't contribute 
to the Kondo screening. Therefore, a dot attached to two metallic leads can be described 
in the framework of the Anderson single-level impurity model. Figure 12.11 illustrates 
schematically one of the spin-fiip co-tunneling processes. The spin-up electron tunnels 
out of the dot, and then it is replaced by the spin-down electron. At low temperature, the 
coherent superposition of all possible co-tunneling processes involving spin-fiip results in 
the screening of the local spin. 



It> 



virtual state 



li> 



i 



4- 



-i- 



Figure 2.1: Spin-fiip cotunneling process in a quantum dot with odd occupation. The left 
and right panels refer to spin- up | t) and spin-down | J.) ground states, which are coupled 
by a cotunneling process. The middle panel corresponds to a high-energy virtual state. 



Let us consider a quantum dot (QD) occupied by one electron with energy Eq in the 
ground state. The dot in tunneling contact with the source and drain leads is described 
by the Anderson Hamiltonian 



Hd + Hlead + Hf 



(2.1) 



Here the first term, Hd, is the Hamiltonian of the isolated dot, 



Hd = y^^Spdld^ + Qn-ini, 



(2.2) 



where a =|, J. is the spin index, Q > is the Coulomb blockade energy, and = d^dfj. 
The second term, Hiead, describes the electrons in the source (s) and drain (d) electrodes. 



Hlead = ^ ^kcl^a^^kaa, a = S,d. 
kcra 

The last term, Ht^ is the tunnehng Hamiltonian, 

Ht = 5^ (^.cLc^. + ff.c.^ 



(2.3) 



(2.4) 



kaa 



where Va {a = s,d) are tunneling matrix elements. Here and below we assume that Va 
are real and positive. It is convenient to perform a canonical transformation 



(2.5) 
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with 



u = — :, V = — , (2.0) 



As a result, only the fermions Cka contribute to tunneling, and the tunneling Hamiltonian 
fl2.4|l takes the form, 

Ht = vY,{ij^ + H.c), (2.7) 



where V = ^V^ + Vj. 

The spectrum of electrons in the leads forms a band with bandwidth 2Dq. In accor- 
dance with the Haldane renormalization group (RG) procedure, the low energy physics 
can be exposed by integrating out the high-energy charge excitations ,54^ . This procedure 
implies the renormalization of the energy level of the dot by mapping the initial energy 
spectrum — -Dq < Cfc < -Dq onto a reduced energy band —-Do + \^^\ < < Dq — \6D\ 
(Fig|23: 

. = .0-^, (2.8) 

where T = npoV"^ is the tunneling rate, po is the density of electron states in the leads, 
which is assumed to be constant. Iterating the renormalization procedure ()2.8|) . one 
obtains the scaling equation, which describes the evolution of the one-electron energy 
state of the dot with reducing the energy scale of the band continuum, 

de r 



dlnD 7t' 



(2.9) 



The processes involving charge scattering to the band edges lead to renormalization of F 
only in higher order in V: 



dlnD 

and hence for T <^ D there is no significant renormalization of F. 
The scaling invariant for equation ()2.9p . 



e{D)-^ln(^]. (2.10) 



TT 



F 



is tuned to satisfy the initial condition e{Do) = Eq. 

The above Haldane RG procedure brings us to the Schrieffer-Wolf (SW) limit [55J D ~ 
|£:(-D)|, where all charge degrees of freedom are quenched for excitation energies within 
the interval —D < < D and scaling terminates. The excited states with two or zero 
electrons are higher in energy by Q+Sq and |£o| {Q+^o, ko| ^ ^), respectively, and are not 
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Figure 2.2: Haldane renormalization group procedure. Reducing the bandwidth by 2|5Z}| 
(panel (a)) results in the renormalization of the energy level of the dot (panel (b)). 

important for the low-energy dynamics of the dot. The effective spin Hamiltonian with the 
two- and zero-electron states frozen out can be obtained by means of the Schrieffer- Wolff 
unitary-transformation ['SF applied to the Hamiltonian ()2.1|) . 



where the operator S is found from the condition 



(2.11) 



(2.12) 



The condition ()2.12|1 means that the effective Hamiltonian ()2.11|) does not contain linear 
in V terms, which allow the variation in the number of electrons in the dot. 

Retaining the terms to order 0{V^) on the right-hand side of Eg. p.llj) . one comes to 
the Kondo Hamiltonian 



J 



Hk = ^kcl^^ka + 7X1 ^kaCk'a + JS ■ S. 



4 

kcT kk'a 

Here S is the spin one-half operator of the dot, 



(2.13) 



[2.U) 



s represents the spin states of the conduction electrons. 



n Y ^k 



:2.i5) 



kk'acr' 



and T is the vector of Pauli matrices. The antiferromagnetic coupling constant is. 



1 



1 



J 

2 \tF - Sq Eq + Q - epj ' 



[2.16) 



Scaling equation for the coupling constant ()2.16|) can be derived by the poor-man's 
scaling method [SG, . The essence of the scaling approach is that the higher energy excita- 
tions can be absorbed as a renormalization of the coupling constant J. To carry out the 
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scaling we divide the conduction band into states, —D + \5D\ < < D — \6D\, which are 
retained, and states within \6D\ of the band edge which are to be ehminated provided the 
effective exchange Hamihonian, JS ■ s, (the last term of the Kondo Hamiltonian ()2.13p ) is 
perturbatively renormalized by changing the coupling constant J ^ J + 6 J. The lowest 
order corrections to J due to virtual scattering of conduction electrons into the band edge 
can be represented by the second order diagrams (Fig. 12. 3|) . Calculating the contribution 
of these diagrams one obtains, 

6J = -poJ'^-^. (2.17) 
Eq. ()2.17p leads to the scaling equation, 

"^"^ -poJ'- (2.18) 



d\n{poD) 

Integrating Eq. ()2.18|1 from the initial band width D and coupling constant J ()2.1(i|l to a 
new band width D and renormalized coupling constant J yields. 



Dexpi = Dexpi ^ . (2.19) 

\ PoJ J V poJj 

Eq. ()2.19|) shows that the solution of the scaling equation ()2.18p is characterized by a 
single parameter which plays the role of a scaling invariant |HZ|. This scaling invariant is 
called Kondo temperature. 



q cr 



^7 ^/,- 
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Figure 2.3: Second order diagrams describing the scattering of a conduction electron from 
the state fcjCTj into an intermediate particle (panel (a)) or hole (panel (b)) state qo" and 
then to a final state fc/O"/. The dashed lines represent the conduction electron, whereas 
the solid lines correspond to the localized spin of the dot. 

The Kondo effect can be observed by measuring the dc current induced by a direct bias 
voltage Vdc applied across the dot. The corresponding differential conductance G'(Vdc) = 
dl/dVdc exhibits a sharp peak at Vdc = 0, which is called zero-bias anomaly (ZBA) 
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[23 1211 inni ini] • At finite value of the bias eVdc ^ the Kondo effect is suppressed 
jlOj . Therefore, the width of the peak of the differential conductance at zero bias is of 
the order of T^- In the weak coupling regime T ^ Tx the Kondo contribution Gk to the 
differential conductance can be calculated by means of Keldysh technique 



37]-2 1 AV'^V^ 
= "l^[ln(r/T^)]2^^' ^^ = ^Jvffv!y- ^^-^^^ 

The Kondo temperature Tk is the only energy scale which controls all low-energy prop- 
erties of the quantum dot. Eq. ()2.2H) shows that the ratio Gk/Gu depends only on the 
dimensionless variable T/Tk jSHj- In the strong coupling regime T <C Tk the spin-flip 
scattering is suppressed, and the system allows an effective Fermi liquid description 
The zero-bias conductance then follows from the Landauer formula, Gk = Gu- 
lf a magnetic fleld is applied to the system, the zero-bias peak splits into two peaks 
at eVdc = ^Ez, where Ez is the Zeeman energy. These peaks are observable even at 
eVdc^Ez ^ Tk [El EHl ED] • However, due to a nonequilibrium-induced decoherence, 
these peaks are wider than Tk, and the value of the conductance at the peaks never reaches 
the unitary limit Gu PII1IS1IS21- In the next section we demonstrate that quantum dots 
with even electron occupation may exhibit a generic Kondo effect at certain value of the 
Zeeman energy Ez ^ Tk- 

2.2 Magnetic Field Induced Kondo Tunneling through 
Evenly Occupied QD 

The Kondo effect discussed in the previous section takes place only if the dot has a nonzero 
spin in the ground state. This is always the case for odd electron occupation A/". When 
J\f is even, the ground state of the spin-degenerate QD is a singlet {S = 0) since all 
single-particle levels are occupied by pairs of electrons with opposite spins. According to 
Hund's rule, the lowest excited state of the dot is a triplet {S = 1) at a distance 6 above 
the ground state. The spacing 6 can be tuned by means of a magnetic fleld. Application 
of a magnetic fleld results in a singlet-triplet transition in the ground state of the dot, 
leading to a Kondo effect. This magnetic field induced Kondo effect occnis both in vertical 
QDs with shell-like structure of electronic states ^1 El QB| and in planar (lateral) QDs 
formed by orbitally non-degenerate electron states ^2]- In the former case the Zeeman 
effect, which lifts the spin degeneracy, is negligibly small in comparison with diamagnetic 
shift because of a small value of the effective (7-factor in semiconductor heterostructures. 
Magnetic fleld affects mostly the orbital states, leading to singlet-triplet level crossing. 
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In the latter case, the Zeeman contribution dominates, and twofold degeneracy of the 
ground state of an isolated dot appears only if the Zeeman energy Ez is equal to the 
single-particle level spacing 5 in the dot. Both types of field-induced Kondo tunneling 
were observed [131 E21 • 

Let us consider a planar QD with even number of electrons, weakly connected to the 
metallic electrodes, and subject to an external magnetic field. In QDs, charge and spin 
excitations are controlled by two energy scales, charging energy Q and single-particle 
level spacing 5 respectively, which typically differ by an order of magnitude: Q ~ ImeV, 
5 ~ O.lmeV [2H1 EHl EUl E] • This separation of energy scales allows one to change the spin 
state of the dot, without changing its charge. In an applied field = S/gfiB l2j, the 
spin-up projection of the triplet becomes degenerate with the singlet ground state (Fig. 
12. 4|) . At this point, the spin-flip transitions shown in Fig. 12.51 become possible, leading 
to a new type of Kondo resonance. It was found ^3] that at a certain magnetic field the 
differential conductance has a peak at zero bias voltage (Fig. 12. 6|) . 




% B 



Figure 2.4: Low-energy states of an evenly occupied QD in magnetic field. The spin- 
up projection |T, 1) of the triplet becomes degenerate with the singlet ground state at 
Be = 5/gnB- 



is> 



virtual state 



H 



4 




|T,1> 



4- 



Figure 2.5: Magnetic field induced Kondo effect in a QD with even occupation. Spin-flip 
transitions connecting the singlet \S) (left panel) and triplet |T, 1) (right panel) states. 
The intermediate high-energy virtual state is shown in the middle panel. 
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Initially, the QD in Fig. 12.51 is treated within an Anderson- type model with bare level 
operators dai-, energies Si and tunneling matrix elements Vi with i = 1, 2 (we consider the 
symmetric case Vis = Vid = V)- Next, the isolated dot Hamiltonian is diagonalized in the 
Hilbert space which is a direct sum of two (|A)), one and three (|A)) electron states, using 
Hubbard operators X'^'^ = |7)(7| (7 = A, A) ^HIIM!- The two particle states |A) exhaust 
the lowest part of the spectrum consisting of a singlet 15) and triplet |T, /i) (/x = 1, 0, — 1). 
The corresponding energies are, 

Es = 2e2, ET^ = ei + e2 + gf^Bf^B, (2.22) 

with El — 62 = S, fiB is the Bohr magneton and g ^ 2 is the free electron g-iactor. The 
energies of one and three electron states are of order of the charging energy. Ex ~ Q- 
Finally, tunneling operators in the bare Anderson Hamiltonian are replaced by a product 
of number changing Hubbard operators X'^^ and a combination Cka = 2^^^'^(ckas + Ckad) 
of lead electron operators, {k =momentum, a = spin projection and s, d stand for source 
and drain). 

With these preliminaries, the starting point is a generalized Anderson Hamiltonian 
describing the QD in tunneling contact with the leads, 

Ha=J211 '^^l.Ckaa + ^7^"" + ( E E + H.c\ , (2.23) 

a=s,d ka 7=AA \ AA ka / 

with dispersion of electrons in the leads and V^^ = Vi{X\dai\A) . The two states of the 
dot which become degenerate at = S/gfiB, are 

|S) = 44,10), |T,1) = 44|0). (2.24) 



Since Q ^ S, one- and three-electron states can be integrated out by means of the 
SW transformation [22]. The resulting exchange Hamiltonian has a form of anisotropic 
Kondo Hamiltonian, 

He. = J\\S,s, + ^ [S+s- + S-S+) . (2.25) 

Here the effective exchange constants are 

2(Vi^ + Vj) AV,V2 
J\\ = Q , = (2.26) 

The spherical components of the dot spin operator S are defined via Hubbard operators 
connecting the \S) and |T, 1) states of the dot, 

5, = ^(Xi'i-X^^), S+ = X''', S- = X''\ (2.27) 
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Figure 2.6: Differential conductance dl/dV as a function of bias voltage V at different 
values of magnetic field B [13j. At -Be = 1-18 T, the differential conductance has a peak 
at zero voltage. 



The conduction electron spin operators are determined by Eq. ()2.15|1 . 

The scaling equations for dimensionless exchange constants ji^ = p^Jy (z^ =||, J-) read, 

- (2.28) 



dhid d\n.d 
yielding the Kondo temperature, 



2j 

where D is the effective bandwidth in the SW limit, and 



Tk = Z^expj - A ) , (2.29) 



J\\-jI \n + ^ J\\ - jI , 

In the isotropic limit j|| — j± ^0 one has A 2 and Eq. ()2.29p reduces to the usual 
expression Tk = -D exp(— l/j||). 



2.3 Singlet- Triplet Kondo Effect in Double Quantum 
Dot 

The isolated quantum dots considered above are typical examples of artificial atoms. A 
double valley quantum dot with weak capacitive and/or tunnelling coupling between its 
two wells may be considered as the simplest case of artificial two-atom molecule. Its 
closest natural analog is the hydrogen molecule H2 or the corresponding molecular ions 
H2 for the occupation number TV" = 2, 1, 3, respectively |65] . 
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d I 

Figure 2.7: Double Quantum Dot in side (T-shaped) geometry. 

Let us consider the DQD with two electron occupation in a side geometry (Fig. I2.7|l . 
Each valley is described by the one-electron level Coulomb blockade energy Qa, and 
bare level operators (ig-a {a = l,r for left and right dot, respectively). The right dot is 
assumed to have a smaller radius and, hence, larger capacitive energy than the left dot, 
i.e., Qr'^Qi. The left dot is coupled by tunneling W to the right dot (W<t^Qi^r) and by 
tunneling V to the source (s) and drain (d) leads. The spectrum of an isolated DQD 
consists of the singlet ground state \S), the low-lying triplet exciton |T, /i) (/i = 1,0,1) 
and high-energy charge transfer singlet exciton \Ex), 

\S) = a\s) - V2(3\ex), 

|T,1) = 4<J0), |T,1)=4/,|0), 



where 




(2.30) 



The corresponding energies are, 

Es = ei + er-2WP, = Q + e,,, Ee^ = 2ei + 2Wi3, (2.31) 

where f3 = W/Air < 1 (A;^ = + q - e^). 

The DQD in tunneling contact with the leads can be described by a generalized An- 
derson Hamiltonian, 

Ha=J211 ''^"i.tCk.t + Yl E,X^' + ( E E + H.cA . (2.32) 

b=s,d ka 7=AA \ AA ka / 

Here |A) are the two-electron eigenfunctions ()2.30p . |A) are the one- and three-electron 
eigenstates; X^^ = |A)(A| are number changing dot Hubbard operators; Vj"^ = V {X\di(j\A) . 
The Kondo effect at T > Tx is unravelled by employing a renormalization group (RG) 
procedure [54j in which the energies E^ are renormalized as a result of rescaling high- 
energy charge excitations (see Eqs. ()2.9|) and ()2.10|) ). Our attention, though, is focused 
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on renormalization of EsiEt ()2.31|) . Since the tunnel constants are irrelevant variables 
|35| inij , the scaling equations are 

dE\ Fa 

TT 



(2.33) 



d\nD 

Here 2D is the conduction electron bandwidth, Fa are the tunneling strengths, 

Ft = TTp^V\ Ts = a'Tr, (2.34) 

with a = a/1 — 2/5^ < 1 and po being the density of states at ep. The scaling invariants 
for equations ()2.33|) . 



VT V i A 



(2.35) 



are tuned to satisfy the initial condition E\{Dq) = E^j^\ Equations ()2.33|1 determine 
two scaling trajectories E/^{D) for singlet and triplet states. Note that the level Eex is 
irrelevant, but admixture of the bare exciton \ex) in the singlet state is crucial for the 
inequality of tunneling rates Fy > F5 [33 EE]- As a result, the energy Et{D) decreases 
with D faster than Es{D), so that the trajectory Et{D, Tt) usually intersects Es{D, F5) 
at a certain point D = D^. This level crossing may occur either before or after reaching 
the Schrieffer- Wolff (SW) limit D where Ex{D) ~ D and scaling terminates [H^- When 
the scaling trajectories cross near the SW boundary ^ D the singlet ground state 
becomes degenerate with a triplet one (Fig. I2.8|l . As a result, the Kondo resonance may 
arise in spite of the even electron occupation of the DQD. 




Figure 2.8: Haldane renormalization group procedure. 

The above Haldane RG procedure brings us to the SW hmit [HS], where all charge 
degrees of freedom are quenched. By properly tuning the SW transformation e*"^ the 
effective Hamiltonian H = e^^HAe~^^ is of the s — d type However, unlike the 

conventional case j^S] of doublet spin 1/2 we have here the degenerate singlet and triplet 
states A = {S,T}, and the SW transformation intermixes these states. To order 0(|yp), 
then, the tunneling exchange Hamiltonian reads. 



(2.36) 



kerb 
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Here the (antiferromagnetic) coupling constants are 

= J'^ = aJ^. (2.37) 

— e« 

The conduction electron spin operator s is defined by Eq. (I2.15j) . S is the dot spin one 
operator with projections /i = 1, 0, 1, while R couples singlet \S) with triplet (T/i|. Their 
spherical components are defined via Hubbard operators: 

S+ = V2(X^° + X°"^), S- = {S+)\ 5^ = X^^-X", 

R+ = V2(Xi^ - X^i), R~ = {R^)\ R' = -(X°^ + (2.38) 

The vector operators R and S obey the commutation relations of 04 Lie algebra, 

[Sj, Sk] = idjkmSmy [-^jy -^k] — '^(^jkmSmi Sk\ — if^jkm Rm. (2.39) 

(here j,k,m are Cartesian indices). Besides, S ■ R = 0, and the Casimir operator is 
+ R^ = 3. The operators S and R manifest the 50(4) dynamical symmetry of the 

DQD. This justifies the qualification of such DQD as a spin rotator (35t l36]. 
Scaling equations for and J^^ are, 

with ji = poJ^,j2 = PqJ^^, d = PqD. If 6 = Et{D) — Es{D) is the smallest energy scale, 
the energy spectrum of the DQD is quasi degenerate and the system ()2.4()j) is reduced to 
a single equation for the effective integral j+ = ji + ^2, 

Then the RG flow diagram has an inflnite flxed point, and the solution of Eq. 1)2.411) gives 
the Kondo temperature 

T^o = ^expf-— ^ . (2.42) 



In the general case, the scaling behavior is more complicated. The flow diagram still 
has a flxed point at inflnity, but the Kondo temperature turns out to be a sharp function 
of 5. In the case 5 < 0, |5| > Tj^, the scaling of J^^ terminates at D ~ |(!| [H UHl ESI Ei • 
Then one is left with the familiar physics of an under-screened spin one Kondo model 
inn]- The flxed point is still at inflnite exchange coupling J^, but the Kondo temperature 
becomes a function of 5, 

J-K -LkO 



T - \ lAI / ' ^2.43) 
-LKO V 1^1 / 

where a < 1 is determined by the DQD parameters (see the text after Eq. ()2.34|1 ). The 
symmetry of the DQD in this case is SO (3). 
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2.4 Dynamical Symmetry or Complex Quantum Dots 



In the previous sections an accidental level degeneracy is induced by an external magnetic 
field (Sec. 12 .21) and the dot-lead interaction (Sec. 12. 3|) . In both cases the nano-object 
possesses dynamical symmetry. In this section we present the concept of dynamical sym- 
metry in more details, and particularly, discuss its emergence in Complex Quantum Dots 
(CQD). 

The term Dynamical Symmetry implies the symmetry of eigenvectors of a quantum 
system forming an irreducible representation of a certain Lie group. The main ideas and 
the relevant mathematical tools can be found, e.g., in Refs. [HZj- Here they are formulated 
in a form convenient for our specific purposes without much mathematical rigor. We have 
in mind a quantum system with Hamiltonian H whose eigenstates |A) = \Mfi) form ( for 
a given M) a basis to an irreducible representation of some Lie group Q. The energies Em 
do not depend on the "magnetic" quantum number /z. For definiteness one may think of 
M as an angular momentum and of fi as its projection, so that Q is just SU{2). Now let 
us look for operators which induce transitions between different eigenstates. An economic 
way for identifying them is through the Hubbard operators 

X^^' = |A)(A'|. (2.44) 

It is natural to divide this set of operators into two subsets. The first one contains the 
operators |M/i)(/i'M| while the second one includes operators |M/i)(/x'M'| (M 7^ M') for 
which \Mfi) and \M'fi') belong to a c?Zj^erent representation space of Q. A central question 
at this stage is whether these operators (or rather, certain linear combinations of them) 
form a close algebra. In some particular cases it is possible to form linear combinations 
within each set and obtain two new sets of operators {S} and {R} with the following 
properties: 1) For a given M the operators {S} generate the M irreducible representation 
of Q and commute with H. 2) For a given set Mj the operators {S} and {R} form an 
algebra (the dynamic algebra) and generate a non-compact Lie group A. The reason for 
the adjective dynamic is that, originally, the operators {R} do not appear in the bare 
Hamiltonian H and emerge only when additional interactions (e.g., dot-lead tunneling) 
are present. In the special case Q = SU{2) the operators in {5*} are the vector S of spin 
operators determining the corresponding irreducible representations, while the operators 
in the set {R} can be grouped into a sequence of vector operators describing transitions 
between states belonging to different representations of the SU (2) group. 

Strictly speaking, the group A is not a symmetry group of the Hamiltonian H since 
the operators {R} do not commute with H. Indeed, let us express H in terms of diagonal 
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Hubbard operators, 



AA 



(2.45) 



K=M^l A 



so that 



{Em — Em')X' 



AA' 



(2.46) 



As we have mentioned above, the symmetry group Q of the Hamiltonian H , is generated 
by the operators x^^^'^'^'^^'^^' . Remarkably, however, the dynamics of CQD in contact 
with metalhc leads and/or an external magnetic field leads to renormalization of the 
energies {Em} in such a way that a few levels at the bottom of the spectrum become 
degenerate, Em^ = Em2 = ■■■EMn- Hence, in this low energy subspace, the group A 
which is generated by the operators {S} and {R} is a symmetry group of H referred 
to as the dynamical symmetry group. The symbol R is due to the analogy with the 
Runge-Lenz operator, the hallmark of dynamical symmetry of the Kepler and Coulomb 
problems. (The Coulomb potential possesses accidental degeneracy of states with different 
angular momentum 1. Hence, according to ()2.46|) the Runge-Lenz vector is an integral 
of the motion. In this case one speaks about hidden symmetry of the system.) Below 
we will use the term dynamical symmetry also in cases where the levels are not strictly 
degenerate, but their differences are bounded by a certain energy scale, which is the Kondo 
energy in our special case. In that sense, the symmetry is of course not exact, but rather, 
approximate. 

Using the notions of dynamical symmetry, numerous familiar quantum objects, such as 
hydrogen atom, quantum oscillator in (i-dimensions, quantum rotator, may be described in 
a compact and elegant way. We are interested in a special application of this theory, when 
the symmetry of the quantum system is approximate and its violation may be treated as a 
perturbation. This aspect of dynamical symmetry was first introduced in particle physics 
[US] , where the classification of hadron eigenstates is given in terms of non-compact Lie 
groups. In our case, the rotationally invariant object is an isolated quantum dot, whose 
spin symmetry is violated by electron tunneling between the dot and leads under the 
condition of strong Coulomb blockade. 

The special cases Q = SU{2) and A = SO{n) or SU{n) is realizable in CQD. Let 
us first recall the manner in which the spin vectors appear in the effective low energy 
Hamiltonian of the QD in tunneling contact with metallic leads. When strong Coulomb 
blockade completely suppresses charge fluctuations in QD, only spin degrees of freedom 
are involved in tunneling via the Kondo mechanism |2S1IZZ1- An isolated QD in this regime 
is represented solely by its spin vector S. This is a manifestation of rotational symmetry 
which is of geometrical origin. The exchange interaction JS • s (s is the spin operator 
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of the metallic electrons) induces transitions between states belonging to the same spin 
(and breaks SU(2) invariance). On the other hand, the low energy spectrum of spin 
excitations in CQD is not characterized solely by its spin operator since there are states 
close in energy, which belong to different representation spaces of SU(2). Incidentally, 
these might have either the same spin S (like, e.g, in two different doublets) or a different 
spin (like, e.g., in the case of singlet-triplet or doublet-quartet transitions). The exchange 
interaction must then contain also other operators (the R-operators mentioned above) 
inducing transitions between states belonging to different representations. The interesting 
physics occurs when the operators R„ "approximately" commute with the Hamiltonian 
Hdot of the isolated dot. In accordance with our previous discussion, the R-operators are 
expressible in terms of Hubbard operators and have only non-diagonal matrix elements 
in the basis of the eigenstates of H^ot- The spin algebra is then a subalgebra of a more 
general non-compact Lie algebra formed by the whole set of vector operators {S,Rri}. 
This algebra is characterized by the commutation relations, 

[5*4, Sj] = itijkSk, [Si, Rnj] = itijkRnk) \Rni-,Rnj\ — '^^ijk^k) (2.47) 

with structure constants tij^, f^-j^ (here ijk are Cartesian indices). The R-operators are 
orthogonal to S, 

S ■ R„ = 0. (2.48) 

In the general case, CQDs possess also other symmetry elements (permutations, reflec- 
tions, finite rotations). Then, additional scalar generators Ap arise. These generators also 
may be expressed via the bare Hubbard operators, and their commutation relations with 
R-operators have the form 

\Rni} Rmj\ ^9ij -^pi [Rnii ^fij -^mj i (2.49) 

with structure constants g^J^^ and /j""*^ {n ^ m). The operators obeying the commutation 
relations ()2.47p and ()2.49|) form an o„ algebra. The Casimir operator for this algebra is 

n p 

Various representations of all these operators via basic Hubbard operators will be estab- 
lished in the following chapters, where the properties of specific CQDs are studied. 

Next, we show how the dynamical symmetry of CQD is revealed in the effective spin 
Hamiltonian describing Kondo tunneling. This Hamiltonian is derived from the general- 
ized Anderson Hamiltonian 

Ha = Hdot + Hiead + Htun- (2-51) 
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The three terms on the right hand side are the dot, lead and tunnehng Hamiltonians, 
respectively. In the generic case, a planar CQD is a confined region of a semiconductor 
secluded between drain and source leads, with complicated multivalley structure. The 
CQD contains several valleys numbered by index a. Some of these valleys are connected 
with each other by tunnel channels characterized by coupling constants Waa', and some 
of them are connected with the leads by tunneling. The corresponding tunneling matrix 
elements are Vab {b = s,d stands for source and drain, respectively). The total number 
of electrons in a neutral CQD as well as the partial occupation numbers jVa for the 
separate wells are regulated by Coulomb blockade and gate voltages Vga applied to these 
wells, with Af — J2a-^a- assumed that the capacitive energy for the whole CQD is 
strong enough to suppress charged states with J\f' — Af ±1, which may arise in a process 
of lead-dot tunneling. 

If the inter- well tunnel matrix elements Waa' are larger than the dot-lead ones Vab (or 
if all tunneling strengths are comparable), it is convenient first to diagonalize Hdot and 
then consider Htun as a perturbation. In this case Hdot may be represented as 

/^<iot=5]^A|A)(A|+ Yl ^a|A)(A|. (2.52) 
AeM XeM±i 

Here all intradot interactions are taken into account. The kets |A) = |jV, q) represent 
eigenstates of Hdot in the charge sector J\f and other quantum numbers q, whereas the 
kets I A) = \J\f ±l,p) are eigenstates in the charge sectors ± 1 with quantum numbers 
p. All other charge states are suppressed by Coulomb blockade. Usually, q and p refer to 
spin quantum numbers but sometimes other specifications are required (see below). 
The lead Hamiltonian takes a form 

Hlead — ^ ] ^kaC^aka'^ocka- (2.53) 
k,a,(T 

In the general case, the individual dots composing the CQD are spatially separated, so 
one should envisage the situation when each dot is coupled by its own channel to the lead 
electron states. So, the electrons in the leads are characterized by the index a, which 
specifies the lead (source and drain) and the tunneling channel, as well as by the wave 
vector k and spin projection cr. 

The tunnel Hamiltonian involves electron transfer between the leads and the CQD, 
and thus couples states |A) of the dot with occupation J\f and states |A) of the dot with 
occupation jVi 1. This is best encoded in terms of non-diagonal dot Hubbard operators, 
which intermix the states from different charge sectors 

X^^ = |A)(A|, X^^ = |A)(A|. (2.54) 
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Thus, 



^*=E E {VaaAj^){M+H.C.) + J2 E {yaaAjX){A\+H.C.), (2.55) 

kaaa AeAT+i kaaa AeAT-l 

where \<4t = ^a(AMa<x|A). 

Before turning to calculation of CQD conductance, the relevant energy scales should 
be specified. First, we suppose that the bandwidth of the continuum states in the leads. 
Da, substantially exceeds the tunnel coupling constants, ^ Waa',Va (actually, we 
consider leads made of the same material with Das = Dad = Dq). Second, each well a in 
the CQD is characterized by the excitation energy defined as = Ex{Afa — 1) — E\{Afa), 
i.e., the energy necessary to extract one electron from the well containing Afa electrons 
and move it to the Fermi level of the leads (from now on the Fermi energy is used as the 
reference zero energy level). Note that is tunable by applying the corresponding gate 
voltage Vga- We are mainly interested in situations where the condition 

A, ~ Do, Q„ (2.56) 

is satisfied at least for one well labelled by the index c. Here Qc is a capacitive energy, 
which is predetermined by the radius of the well c. Eventually, this well with the largest 
charging energy is responsible for Kondo-like effects in tunneling, provided the occupation 
number Afc is odd. The third condition assumed in most of our models is a weak enough 
Coulomb blockade in all other wells except that with a = c, i.e., Qa ^ Qc- Finally, we 
demand that 

baa = ^ « 1, (2-57) 
for those wells, which are coupled with metallic leads, and 

/3a = ^ « 1. (2.58) 



-'ac 



Here Eac are the charge transfer energies for electron tunneling from the c-well to other 
wells in the CQD. 

The interdot coupling under Coulomb blockade in each well generates indirect ex- 
change interactions between electrons occupying different wells. Diagonalizing the dot 
Hamiltonian for a given Af = Yla-^a^ easily finds that the low-lying spin spectrum 
in the charge sectors with even occupation A/" consists of singlet /triplet pairs (spin S = 
or 1, respectively). In charge sectors with odd A/" the manifold of spin states consists of 
doublets and quartets (spin S'=l/2 and 3/2, respectively). 

The resonance Kondo tunneling is observed as a temperature dependent zero bias 
anomaly in tunnel conductance [2H1 122] • According to existing theoretical understanding. 
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the quasielastic cotunneling accompanied by the spin flip transitions in a quantum dot is 
responsible for this anomaly. To describe the cotunneling through a neutral CQD with 
given A/", one should integrate out transitions involving high-energy states from charge 
sectors with M' = A/'± 1. In the weak coupling regime at T > this procedure is done 
by means of perturbation theory which can be employed in a compact form within the 
renormalization group (RG) approach formulated in Refs. [511 EEl- 

As a result of the RG iteration procedure, the energy levels in the Hamiltonian 
()2.52|) are renormalized and indirect exchange interactions between the CQD and the 
leads arise. The RG procedure is equivalent to summation of the perturbation series at 
T > Tki where is the Kondo energy characterizing the crossover from a perturbative 
weak coupling limit to a non-perturbative strong coupling regime. The leading logarithmic 
approximation of perturbation theory corresponds to a single-loop approximation of RG 
theory. Within this accuracy the tunnel constants W and V are not renormalized, as well 
as the charge transfer energy ()2.56|) . Reduction of the energy scale from the initial value 
Dq to a lower scale ~ T results in renormalization of the energy levels i?A E\{Do/T) 
and generates an indirect exchange interaction between the dot and the leads with an 
(antiferromagnetic) exchange constant J. 

The rotational symmetry of a simple quantum dot is broken by the spin-dependent 
interaction with the leads, which arises in second order in the tunneling amplitude Va- 
in complete analogy, the dynamical symmetry of a composite quantum dot is exposed 
(broken) as encoded in the effective exchange Hamiltonian. In a generic case, there 
are, in fact, several exchange constants arranged within an exchange matrix J which 
is non-diagonal both in dot and lead quantum numbers. The corresponding exchange 
Hamiltonian is responsible for spin-flip assisted cotunneling through the CQD as well as 
for singlet-triplet transitions. 

The precise manner in which these statements are quantified will now be explained. 
After completing the RG procedure, one arrives at an effective (or renormalized) Hamil- 
tonian H in a reduced energy scale D, 



H — Hdot + Hiead + Hcotuni 



(2.59) 



where the effective dot Hamiltonian ()2.52p is reduced to 




(2.60) 



written in terms of diagonal Hubbard operators. 



X 



AA 



|A)(A|. 



(2.61) 
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At this stage, the manifold {A} G M contains only the renormalized low-energy states 
within the energy interval comparable with Tk (to be defined below). Some of these states 
may be quasi degenerate, with energy differences \E{^ — Ex'\ < Tk- However, Tk itself 
is a function of these energy distances (see, e.g., [121 IMl HO]), and all the levels, which 
influence Tk, should be retained in ()2.(j()j) . 

The effective cotunneling Hamiltonian acquires the form 

Hcot = Yl ( ^0""'S ■ S""' + J2 ^n"'Rn " s""' J . (2.62) 
aa' \ n / 

Here S is the spin operator of CQD in its ground state, the operators s""' represent the 
spin states of lead electrons, 

kk' crcr' 

where f is the vector of Pauli matrices. In the conventional Kondo effect the logarithmic 
divergent processes develop due to spin reversals given by the first term containing the 
operator S. In CQD possessing dynamical symmetry, all R- vectors are involved in Kondo 
tunneling. In the following chapters we will show how these additional processes are 
manifested in resonance Kondo tunneling through CQD. Note that the elements of the 
matrix J are also subject to temperature dependent renormalization J^"' J^"' (Dq/T). 

The cotunneling Hamiltonian ()2.62j) is the natural generalization of the conventional 
Kondo Hamiltonian JS ■ s for CQDs possessing dynamical symmetries. In many cases 
there are several dot spin 1 operators depending on which pair of electrons is "active". 
In this pair, one electron sits in well c and the other one sits in some well a. The other 
M —2 electrons are paired in singlet states. This scenario applies if M is even. The spin 1 
operator for the active pair is denoted as Sq. (In some sense, the need to specify which pair 
couples to S' = 1 while all other pairs are coupled to S' = is the analog of the seniority 
scheme in atomic and nuclear physics (see, e.g., plj)). The cotunneling Hamiltonian for 
CQD contains exchange terms J^"'Sa-s""'. Then, instead of a single exchange term (first 
term on the RHS of Eq. ()2.62|l ). one has a sum J2"'^a ■ s°"'. Additional symmetry 
elements (finite rotations and refiections) turn the cotunneling Hamiltonian even more 
complicated. In the following chapters we will consider several examples of such CQDs. 
It is seen from ()2.62j) . that in the generic case, both spin and R- vectors may be the 
sources of anomalous Kondo resonances. The contribution of these vectors depends on 
the hierarchy of the energy states in the manifold. In principle, it may happen that the 
main contribution to the Kondo tunneling is given not by the spin of the dot, but by one 
of the R-vectors. 
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Thus, we arrive at the conclusion that the regular procedure of reducing the full 
Hamiltonian of a quantum dot in junctions with metallic leads to an effective Hamiltonian 
describing only spin degrees of freedom of this system reveals a rich dynamical symmetry 
of CQD. Strictly speaking, only an isolated QD with A/" = 1 is fully described by its 
spin 1/2 operator obeying SU{2) symmetry without dynamical degrees of freedom. Yet 
even the doubly occupied dot with M = 2 possesses the dynamical symmetry of a spin 
rotator because its spin spectrum consists of a singlet ground state (S) and a triplet 
excitation (T). Therefore, an R- vector describing S/T transitions may be introduced, and 
the Kondo tunneling through a dot of this kind may involve spin excitation under definite 
physical conditions, e.g., in an external magnetic field A two-electron quantum dot 
under Coulomb blockade constitutes apparently the simplest non-trivial example of a 
nano-object with dynamical symmetry of a spin rotator possessing an SO{A) symmetry. 

Dynamical symmetries SO{n) of CQDs are described by non-compact semi-simple 
algebras [72]. This non-compactness implies that the corresponding algebra o„ may be 
presented as a direct sum of subalgebras, e.g., 04 = 03 © 03. Therefore, the dynamical 
symmetry group may be represented as a direct product of two groups of lower rank. In 
case of spin rotator the product is 5*0(4) = SU{2)®SU (2). Generators of these subgroups 
may be constructed from those of the original group. The 5*0(4) group possesses a single 
R-operator R, and the direct product is realized by means of the transformation 

K=^, N=^. (2.64) 

2 ' 2 ^ ' 

Both vectors K and N generate SU{2) symmetry and may be treated as fictitious S=l/2 
spins [n^. In some situations these vectors are real spins localized in different valleys of 
CQD. In particular, the transformation ()2.64|1 maps a single site Kondo problem for a 
DQD possessing 50(4) symmetry to a two-site Kondo problem for spin 1/2 centers with 
an SU (2) symmetry (see discussion in Refs. P^IIH^ ). For groups of higher dimensionality 
{n > 4) one can use many different ways of factorization, which may be represented by 
means of different Young tableaux (see Appendix D). 

Even in the case n = 4, the transformation ()2.(i4|l is not the only possible "two-spin" 
representation. An alternative representation is realized in an external magnetic field jSHl • 
When the ground state of S/T manifold is a singlet (the energy 6 = Et — Es > 0), the 
Zeeman splitting energy of a triplet in an external magnetic field may exactly compensate 
the exchange splitting 6. This accidental degeneracy is described by the pseudospin 1/2 
formed by the singlet and the up projection of spin 1 triplet. Two other projections of the 
triplet form the second pseudospin 1/2. The Kondo effect induced by external magnetic 
field observed in several nano-objects was the first experimental manifestation of 
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dynamical symmetry in quantum dots. 

In conclusion, we outlined in this section the novel features which appear in effective 
Kondo Hamiltonians due to the dynamical symmetry of CQD exhibiting Kondo tunneling. 
In the following chapters we will see how the additional terms in the Hamiltonian ()2.62j) 
influence the properties of Kondo resonance in various structures of CQDs. 
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Chapter 3 

Trimer in Parallel Geometry 



This Chapter is devoted to a systematic exposure of the Kondo physics in evenly occupied 
artificial trimer, i.e., triple quantum dot (TQD) in parallel geometry. We are interested 
in answering questions pertaining to the nature of the underlying symmetry of the trimer 
Hamiltonian and the algebra of operators appearing in the exchange Hamiltonian. The 
energy spectrum of the isolated TQD is discussed in Section In Section \3.^A the 

renormalization group equations are derived and various cases of accidental degeneracy 
arising due to dot-lead interaction are discussed. It is shown that the TQD manifests 
SO{n) dynamical symmetry in Kondo tunneling regime. The effective spin Hamiltonians 
are written down and the corresponding On algebras are constracted for the P x 5*0(4) x 
5*0(4), 5*0(5) and 5*0(7) dynamical symmetries in Subsections \3.2.l\ \S.2.Si and \S.2.^ 
respectively. The scaling equations are derived and the Kondo temperatures are calculated 
for the cases of P x 5*0(4) x 5*0(4) and 5*0(5) symmetries. The results are summarized 
in the Conclusions. 

3.1 Energy Spectrum 

Double quantum dot with occupation Af = 2 discussed in Sec. 12 .HI is the analog of a 
hydrogen molecule in the Heitler-London limit jSHlEni! and its 5*0(4) symmetry reflects 
the spin properties of ortho/parahydrogen. A much richer artificial object is a triple 
quantum dot (TQD), which can be considered as an analog of a linear molecule RH2. 
The central (c) dot is assumed to have a smaller radius (and, hence, larger capacitive 
energy Q^) than the left (/) and right (r) dots, i.e., Qc ^ Qi^r- Fig. 13.11 illustrates this 
configuration in a parallel geometry, where the " left-right" (/ — r) reflection plane of the 
TQD is perpendicular to the " source- drain" (s — d) reflection plane of metallic electrodes. 
To regulate the occupation of TQD as a whole and its constituents in particular, there 
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Figure 3.1: Triple quantum dot in parallel geometry and energy levels of each dot Ea = 
Ca — Vga (bare energy minus gate voltage). 

is a couple of gates Vgi,Vgr applied to the l,r dots. The energy levels of single- and two- 
electron states in each one of the three constituent dots are shown in the lower panel of 
Fig. 13.11 Here the gate voltages Vgi^r are applied in such a way that the one-electron level 
ec of a c-dot is essentially deeper than those of the /, r-dots, so that the condition ()2.56p 
is satisfied for the c dot, whereas the inequalities ()2.57j) and ()2.58j) are satisfied for the 
"active" / and r dots. Tunneling between the side dots l,r and the central dot c with 
amplitudes Wi^r determines the low energy spin spectrum of the isolated TQD once its 
occupation A/" is given. This system enables the exposure of much richer possibilities for 
additional degeneracy relative to the DQD setup mentioned above due to the presence of 
two channels {l,r). 

The full diagonalization procedure of the Hamiltonian Hdot for the TQD is presented 
in Appendix El When the condition ()2.58j) is valid, the low-energy manifold for A/" = 4 is 
composed of two singlets \Si), \Sr), two triplets \Ta) = \fia) = l,^,fJ'a = la, 0^, !«) and 
a charge transfer singlet exciton \Ex) with an electron removed from the c-well to the 
"outer" wells. Within first order in <^ 1 the corresponding energies are, 

Es^ = ec + ea + 2ea + Qa-2Wa/3a, 

Et, = e, + ea + 2ea + Q-a, (3.1) 
Ee, = 2ei+2er + Qi + Qr + 2Wipi + 2Wrf3r, 

where the charge transfer energies in Eq. ()2.58|) (for determining /?«) are Eac = Qa + ^a — ^c', 
the notation a = l,r and a = r, / is used ubiquitously hereafter. 
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The completely symmetric configuration, ei = = e, Qi = Qr = Q, Wi = W,. = 
should be considered separately. In this case the singlet states form even and odd 

combinations in close analogy with the molecular states in axisymmetric molecules. 

The odd state S_ and two triplet states are degenerate: 

Es+=ec + ^e + Q- 4Wf3, 

Es- = Et, = + 3e + g, (3.2) 
Ee. = 4e + 2Q + AW 13. 

Consideration of these two examples provide us with an opportunity to investigate the 
dynamical symmetry of CQD. 

3.2 Derivation and Solution of Scaling Equations 

We consider the case of TQD with even occupation AT = 4 discussed in Refs. [731 Ull • 
This configuration is a direct generalization of an asymmetric spin rotator, i.e., the double 
quantum dot in a side-bound geometry |^. Compared with the asymmetric DQD, this 
composite dot possesses one more symmetry element, i.e., the / — r permutation, which, 
as will be seen below, enriches the dynamical properties of CQD. 

Following a glance at the energy level scheme ()3.1|1 , one is tempted to conclude outright 
that for finite W, the ground state of this TQD configuration is a singlet and consequently 
there is no room for the Kondo effect to take place. A more attentive study of the tunneling 
problem, however, shows that tunneling between the TQD and the leads opens the way 
for a rich Kondo physics accompanied by numerous dynamical symmetries. 

Indeed, inspecting the expressions for the energy levels, one notices that the singlet 
states Es^ are modified due to inter-well tunneling, whereas the triplet states Et^ are 
left intact. This difference is due to the admixture of the singlet states with the charge 
transfer singlet exciton (see Appendix^. As was mentioned in the previous chapter, the 
Kondo cotunneling in the perturbative weak coupling regime at T, e > Tk is excellently 
described within RG formalism [JH EH] • According to general prescriptions of this the- 
ory, the renormalizable parameters of the effective low-energy Hamiltonian in a one-loop 
approximation are the energy levels E^ and the effective indirect exchange vertices J^A'- 

3.2.1 P X 50(4) X 50(4) Symmetry 

To apply the RG procedure to the Kondo tunneling through TQD, let us first specify 
the terms -ff/ead and Htun in the Anderson Hamiltonian ()2.5H) . The most interesting 
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for us are situations where the accidental degeneracy of spin states is reahzed. So we 
consider geometries where the device as a whole possesses either complete or slightly 
violated / — r axial symmetry. Then the quantum number a in -ff/ead (I2.53|) contains 
the lead index (s, c?) and the channel index (/,r). The two tunneling channels are not 
independent because of weak interchannel hybridization in the leads. This hybridization 
is characterized by a constant t/,. <^ Dq, which is small first due to the angular symmetry, 
and second due to significant spatial separation between the two channels. The wave 
vector k is assumed to remain a good quantum number. Then, having in mind that in 
our model euas = ^kad = ^ka, the generalized Hamiltonian ()2.53|1 acquires the form 

Hlead = (^^ka^abka + t IrClbka^abkcr^ ■ (3.3) 
fccr b=s,da=l,r 

The tunneling Hamiltonian ()2.55j) is written as 

Htun = J2I1 T.^^abtclbk.X''' + H.C.). (3.4) 

AA ka ab 

We assume below Vas = Vad = K (see Fig 13. Ij) . 
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Figure 3.2: RG diagrams for the energy levels i?A (o) and the effective exchange vertices 
•^AA' (^) (s6e text for further explanations). 

The iteration processes, which characterize the two-step RG procedure contributing 
to these parameters are illustrated in Fig. 13.21 The intermediate states in these diagrams 
are the high-energy states \q) near the ultraviolet cut-off energy D of the band continuum 
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in the leads (dashed hnes) and the states |A) G A/" — 1 from adjacent charge sectors, which 
are admixed with the low-energy states |A) G A/" by the tunneling Hamiltonian Ht ()2.55|) 
(full lines). For the sake of simplicity we confine ourselves with three-electron states in 
the charge sector. 

In the upper panel, the diagrams contributing to the renormalization of Hdot are 
shown. In comparison with the original theory [54 , this procedure not only results in 
renormalization of the energy levels but also an additional hybridization of the states 
|Aa) via channel mixing terms in the Hamiltonian ()3.4|) . Due to the condition ()2.56|) . the 
central dot c remains "passive" throughout the RG procedure. 

The mathematical realization of the diagrams displayed in Fig. 3.2a is encoded in the 
scahng equations for the energy levels £'a, 

vrd^A Fa 

lilT^Eun^- (3-5) 

Here i?AA = E/^ — Ex, Fa are the tunnel coupling constants which are different for different 
A, 

Tt. = vrpo(K' + 2Vi), Ts^ = alTT^. (3.6) 



Here aa = a/1 — 2/3^, and po is the density of electron states in the leads, which is 
supposed to be energy independent. These scaling equations should be solved at some 
initial conditions 

E^{Do) = E^^\ (3.7) 

where the index (0) marks the bare values of the model parameters entering the Hamil- 
tonian Ha dzni)- 

Besides, the diagram of Fig. 3.2a generates a new vertex M^^' , where the states A, A' 
are either two singlets 5*^,5'^ or two triplets Ti,Tr. The third order Haldane iteration 
procedure results in a scaling equation, 

dMir 7 



dD D 



2 



(3.^ 



with an initial condition Mir{Do) = and a flow rate 7 = poViVrtir- After performing the 
Haldane procedure we formally come to the scaled dot Hamiltonian 

H^t = Yl Ea^X^"""" + MirX'-'^'-' (3.9) 

Aa AaAa 

with the parameters E^^ and Mir depending on the running variable D. 

Due to the above mentioned dependence of tunneling rates on the index A, namely the 
possibility of F^ > F5 and Ts_ > ^s+, the scaling trajectories E\{D) may cross at some 
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value of the monotonically decreasing energy parameter D. The nature of level crossing 
is predetermined by the initial conditions ()3.7|1 and the ratios between the tunneling rates 
Fa. As long as the inequality \E\x\ ^ D is effective and all levels are non-degenerate, 
the scaling equations ()3.5|) may be approximated by Eqs. ()2.33|) . The scaling trajectories 
are determined by the scaling invariants ()2.35j) for equations fl3.5|l . tuned to satisfy the 
initial conditions. With decreasing energy scale D these trajectories flatten and become 
D-independent in the so called Schrieffer- Wolff (SW) limit, which is reached when the ex- 
citation energies become comparable with D. The corresponding effective bandwidth 
is denoted as D (we suppose, for the sake of simplicity, that Aa < Qa, so that only the 
states I A) with A/"' = A/" — 1 are relevant). The simultaneous evolution of interchannel 
hybridization parameter is described by the solution of scaling equation ()3.8|1 . 



If this remarkable level crossing occurs at D > Z), we arrive at the situation where 
adding an indirect exchange interaction between the TQD and the leads changes the mag- 
netic state of the TQD from singlet to triplet. Those states E/^, which remain close enough 
to the new ground state are involved in Kondo tunneling. As a result, the TQD acquires a 
rich dynamical symmetry structure instead of the trivial symmetry of spin singlet prede- 
termined by the initial energy level scheme 1)3.11) . Appearance of the enhancement of the 
hybridization parameter Mir ()3.10|) does not radically influence the general picture, pro- 
vided the flow trajectories cross far from the SW line , due to a very small hybridization 
7 ^ Fa ^ Z). However, we are interested just in cases when the accidental degeneracy 
occurs at the SW line. Various possibilities of this degeneracy are considered below. 

The flow diagrams leading to a non-trivial dynamical symmetry of TQD with A/" = 4 
are presented in Figs. 13.31 13.51 13.61 The horizontal axis on these diagrams corresponds to 
the dimensionless energy scale D/Dq for lead electrons, where the vertical axes represent 
the energy levels Ej>^{D). The dashed line E = —D establishes the SW boundary for these 
levels. 

Before turning to highly degenerate situations, where the system possesses specific 
SO{n) symmetry, it is instructive to consider the general case, where all flow trajectories 
Ej>^{D) are involved in Kondo tunneling in the SW limit. This happens when the whole 
octet of spin singlets and triplets forming the manifold 1)3.1)) remains within the energy 
interval ~ in the SW limit. The level repulsion effect does not prevent the formation 
of such multiplet, provided tir is small enough and the inequality 




(3.10) 



Mir{D) < Tk 



(3.11) 
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is valid. At this stage, the SW procedure for constructing the effective spin Hamiltonian in 
the subspace Ms = {Ti, Si,Tr, Sr} should be applied. This procedure excludes the charged 
states generated by Hf to second order in perturbation theory (see, e.g., |F7j). 

The effective cotunneling Hamiltonian can be derived using Schrieffer-Wolf procedure 
(see Appendix EI). To simplify the SW transformation, one should first rationalize 
the tunneling matrix V in the Hamiltonian ()3.4|) . This 4x4 matrix is diagonalized in the 
s — d,l — r space by means of the transformation to even/odd combinations of lead electron 
fc-states and similar symmetric/antisymmetric combinations of l,r electrons in the dots. 
The form of this transformation for symmetric TQD can be found in Appendix El Like 
in the case of conventional QD [211, this transformation eliminates the odd combination 
of s — d electron wave functions from tunneling Hamiltonian. 

It should be emphasized that this transformation does not exclude the odd component 
from Htun in case of TQD in a series geometry . The same is valid for the Hamiltonians 
()3.3|) . ()3.4|) with tir = 0: in this case the rotation in s — d space conformally maps the 
Hamiltonian Ha ()2.5H1 for TQD in parallel geometry onto that for TQD in series. Both 
these cases will be considered in Chapter H] 

Unlike the case of DQD studied in Refs. [HSl IMl, where the spin operators are the 
total spin S and a single R-operator, describing S/T transitions, the TQD is represented 
by several spin operators corresponding to different Young tableaux (see Appendix D). 
To order 0{\V\'^), then. 



Aa AaA(j ka b=s,da=l,r 

+ 5^ JJS„ ■ s, + J^.P ^ S„ ■ Sa„ + ^ jf^R, ■ + J,,, ^ R„ ■ s„,. (3.12) 



Here we recall that Ej>^^ = E\^{D), Mij. = Mir{D), and the effective exchange constants 
are 



a=l,r 



a=l,r 



a=l,r 



a=l,r 





J: 



-ST 





(3.13) 



a 



The vector operators S^, Rq, R^ and the permutation operator P manifest the dynamical 
symmetry of TQD in a subspace Mg- The permutation operator 




(3.14) 



a=l,r 



/i=l,0,l 



commutes with Si + Sr and R/ + R,.. 
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The spherical components of these vectors are defined via Hubbard operators connect- 
ing different states of the octet, 

R- = {Rt)\ K = -(X°«^" + X^'-O"), (3.15) 
R- = iRt)\ K = -{a-aX'^'-- + 

In addition to the spin operator ()2.(')3|1 for conduction electrons, new spin operators are 
required, 

Snn. = - > > cl,^^f„rT'Cnh'rT'. (3.16) 



Rn 



v/2(X^°°" +X°"^"), 
V2{aaX^''^-^ - a^X^"^^) 



2 '^Lcr^o-cr'Cafc'cr'- 
kk' era' 



An extra symmetry element (1-r permutation) results in more complicated algebra which 
involves new R-operator R and the permutation operator P interchanging / and r com- 
ponents of TQD. 

One can derive from the generic Hamiltonian ()3.12|1 more symmetric effective Hamilto- 
nians describing partly degenerate configurations illustrated by the fiow diagrams of Figs. 
13 .31 13. 5| I3.fil These are the cases when the level crossing occurs in the nearest vicinity of 
the SW line in the fiow diagram. It is important to distinguish between the cases of generic 
and accidental symmetry. In the former case the device possesses intrinsic / — r and s — d 
symmetry, i.e., the left and right dots are identical, the corresponding tunnel parameters 
are equal, and left and right leads also mirror each other, namely, eki = ekr = Cfc- In the 
latter case the gate voltages violate / — r symmetry, e.g., they make e; 7^ e^, VJ 7^ Vr, etc. 
The level degeneracy is achieved due to competition between the / — r interdot tunneling 
and the lead-dot tunneling without changing the symmetry of the Hamiltonian. 




Figure 3.3: Scaling trajectories for P x S'0(4) x S0{4) symmetry in the SW regime. Inset: 
Zoomed in avoided level crossing pattern near the SW line. 



The basic spin Hamiltonian ()3.12|1 acquires a more compact form, when a TQD pos- 
sesses generic or accidental degeneracy. In these cases the operators ()3.15p form close 
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algebras, which predetermine the dynamical symmetry of Kondo tunneling. We start 
the discussion of the pertinent SO (n) symmetries with the most degenerate configuration 
(Fig. 13. 3|) . where the TQD possesses generic / — r axial symmetry, i.e., the left and right 
dots are completely equivalent. Then the energy spectrum of an isolated TQD is given 
by Eqs. ()3.2j) . The four-electron wave functions are calculated in Appendix (^. Such 
TQD is a straightforward generalization of the so called T-shaped DQD introduced in 
Refs. P3 EEl IZni EZl EH] • It is clear, that attachment of a third dot simply adds one more 
element to the symmetry group SO (4), namely the / — r permutation P, which is parity 
sensitive. 

To reduce the Hamiltonian ()3.12|) into a more symmetric form, we rewrite the Hubbard 
operators in terms of new eigenstates ^a, recalculated with taking account of the generic 
degeneracy ()3.2j) and / — r mixing M^^. In assuming that the latter coupling parameter 
is the smallest one, it results in insignificant additional remormalization ~ ^\Mir\'^/{6 + 
Q — Ec) of the states Es+ and Eex- Besides, it intermixes the triplet states and changes 
their nomenclature from left/right to even/odd. The corresponding energy levels are 

ET±{D)=ET^TMir. (3.17) 

The flow trajectories for two pairs of states (T+, T_) and {S+, S-) diverge slowly with 
decreasing D. If this divergence is negligible in the scale of Tk, then three nearly coincident 
trajectories Et±,Es_ cross the fourth trajectory Es+ at some point, since the inequality 
< TTi = Ts- with Tri = SnpoV^, Ts^ = aFri is vahd (a = ^/l - 4/52 < 1). If this 
level crossing happens near the SW line, we arrive at a case of complete degeneracy of the 
renormalized spectrum, and the whole octet Ms is involved in the dynamical symmetry 
(Fig. 13. 3p . The fine structure of the flow diagram in the region of avoided level crossing 
is shown in the inset. 

Since the tunneling occurs in even and odd channels independently, the parity is 
conserved also in indirect SW exchange. As a result, the effective spin Hamiltonian ()3.12j) 
acquires the form 

ka r]=g,u »?=9,m r)=9,u 

-\- J2 ^ ^ Syfj ■ S^fj + ^ ^ {^27] -f^r;ij "I" ^27] '^r/^) ' ^VV' (3.18) 

r]=g,u r)=g,u 

Here ekg = ^k — tir, ^ku = ^k + Ur and the lead operators c^fco- ij] = 9, u) are defined in 
Appendix IHl The operators S^, are defined analogously to Sq, Ra in Eq. ()3.15|) . and 
the vector operators S^^^, R^,^'', R^l'* are defined as: 

^r]fj ~ X^^^fj^ R-ry?^ ~t~ ~ ^^^^R.f^. (3.19) 
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The spherical components of the operators lV„ij and lV„ij are given by 



= y2X^^^^ 4f = (4f+)t, 4f = -X°^^^ (3.20) 



The spin operators for the electrons in the leads are introduced by the obvious relations 

fefc' ctct' fefc' era' 

^gu = 2 ^ ^ C^^^Tg-g-zC^fc/Q-', S„g = (Sp„)^, (3-21) 

instead of ()2.63|) . Now the operator algebra is given by the closed system of commutation 
relations which is a generalization of the 04 algebra, 

The operators are orthogonal to R^, and the Casimir operators in this case are /C^ = 
S,^ + R,^ = 3. This justifies the qualification of such TQD as a double spin rotator which 
is obtained from the spin rotator considered in Refs. [H^l IHB] by a mirror reflection. The 
symmetry of such TQD is P x 50(4) x 50(4). 

Four additional vertices appear in the effective spin Hamiltonian ()3.18j) at the second 
stage of Haldane- Anderson scaling procedure jSHI- As a result, the exchange part of the 
Hamiltonian ()3.18p takes the form 



Hcot — ^ ^ Jlrj^ri ' ^rj + ^ ] Jlrj Rr; ' ^r] ~^ ^2 ^ ] ^VV ' ^VV 

r)=g,u r) r) 

+ ^^(</a7 -R-r/jJ + J2fi '^\^f|) ' ^nV + "^Sy^V ' + ^iq ' ^f}' (3.23) 

The coupling constants in the Hamiltonian ()3.23p are subject to renormalization. Their 
values at D = D are taken as initial conditions 

Jl,{D) = Jl{D) = J^:{D) = Jl^{D) = ^, (3.24) 

ep — £ 

JUD) = J.^JiD) = 0, Jf (D) = a JJ(D) (z = 1, 2, 3) 
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for solving the scaling equations. These can be written in the following form: 



dj 



Ir, 



dlnd 



(aT\2 (jST\2^ 
■T\2 I o( aT I ( aST\2 , \J2 ) , \J2ri) 



+ 



dlnd 

dlnd 
dlnd 

dji-q 

dlnd 

djil 

dlnd 



9 ^-^2 Ulri + Jsr]) + j2r] Ul-q + jfr? )} 



r]=g,u 



- I SjiVfJ + 2(-l)''m,,jfJ + jIjI, 



UP' 



+ 



^^ST,2 



/2fj 



(3.25) 



where jir^ = poJir^ {i = 1,2,3), d = PqD and m^^ = PoMi^. It should be noted that the 
terms proportional to uiir arise in Eqs. ()3.25p since the dot Hamiltonian (the first term in 
Eq. ()3.18p ) is not proportional to the unit matrix, and thus it does not commute with the 
exchange terms ()3.23|1 . As can be seen from Eq. ()3.17|l . the deviation from the unit matrix 
is proportional to Mi^. 

Solution of Eqs. ()3.25p yields the Kondo temperature 

1 



-KO 



D\l- 



2m 



Ir 



(v^+i)(3j?;+jif 

The limiting value of this relation for independent /, r channels is 

4 



(3.26) 



lim Tko = D exp 



iV3 + mjf^ + jfp 



(3.27) 



Here and below the coupling constants ji{D) in all equations for Tx are taken at D = D. 
We see that avoided crossing effect in the case of slightly violated I — r symmetry of TQD 
turns the Kondo temperature to be a function of the level splitting ()3.17|) . A similar 
situation has been noticed in previous studies of DQD [231 EHl and planar QD with even 
occupation [JS], where Tk turned out to be a monotonically decreasing function of S/T 
splitting energy 6 = Eg — Ex with a maximum at 5 = 0. Now the Kondo temperature is a 
function of two parameters, TK{Mir,6). Looking at Fig l3.3l (which corresponds to 5 = 0) 
we notice that for large enough Mir, when the inequality ()3.11|) is violated, Mir ^ T^, 
the symmetry of TQD is reduced to 5*0(4) symmetry of S/T manifold with the Kondo 
temperature 



T, 



Kl 



D exp 



(3.28) 
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Tko 4 

Figure 3.4: Variation of Tk with the parameters 6 and Mir (see text for further details). 

Additional S/T splitting induced by the gate voltage (e^ ^ e^) results in further decrease 
of Tk as a function of 6. The asymptotic form of the function Tk{S) is 



Tk I Tki 



Tki \ S 



(3.29) 



where a = -^/l — 4/5^ < 1. In the limit of 6 ^ D the singlet state should be excluded 
from the manifold, and the symmetry of the TQD with spin one in this case is S0{3). 
The general shape of TK{Mir,6) surface is presented in Fig. 13.41 Thus the Kondo effect 
for the TQD with mirror symmetry is characterized by the stable infinite fixed point 
characteristic for the underscreened spin one dot, similar to that for DQD ISSUSE!- 

3.2.2 S0{5) Symmetry 

Now we turn to asymmetric configurations where Ei^ ^ Ere, Tt^ 7^ Tt;- In this case the 
system loses the l — r symmetry, and it is more convenient to return to the initial variables 
used in the generic Hamiltonian ()3.12|) . 

When the Haldane renormalization results in an accidental degeneracy of two singlets 
and one triplet, Esi ~ -E^, ~ Es^ < Et,. (Fig. 13. 5|) . the TQD acquires an 5*0(5) symmetry 
of a manifold {T;, Si, Sr}- In this case the SW Hamiltonian ()3.12|1 transforms into 

H = Yl E^X^'' + Mir{X'''^+X'^^^)+J2J2('^-''ak.Caka + tlr^^^^ 

A=Ti,Si,Sr kcr a=l,r 

+ JiSz-s, + J2RrSi + J3(Ri-s,/ + R2-Si,), (3.30) 

where Ji = Jf' , J2 = J^^ and J3 = arJh-- The spherical components of the vector 
operators Ri and R2 are given by the following expressions, 

R+ = (/^^)^ R^ = {Rt)\ R2. = Rl. (3.31) 
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Figure 3.5: Scaling trajectories resulting in an 5*0(5) symmetry in the SW regime. Inset: 
Zoomed in avoided level crossing pattern near the SW line. 

The group generators of the 05 algebra are the /-vectors Si, Hi from ()3.15|1 and the oper- 
ators intermixing /- and r-states, namely the vector R = Ri + R2, 

^+ = V2(X^''^'- - X^'-^O, R~ = (^^)^ = + X^'-O'), (3.32) 

and a scalar A interchanging l,r variables of degenerate singlets 

A = i{X^^^' - X^'^-) . (3.33) 

The commutation relations ()2.47|) . ()2.49|) in this particular case acquire the form 

[Slj,Sik] = i^jkmSlmi [RljiRlk\ — i^jkmSlmi {Rlj-,Slk\ = i^jkmRlmi 

\Rj,Rk\ = i^jkmSlmy [Rj^Sn^] = iejkmRm: [Rlj : Rk] = ^^jkA, 

[Rj,A] = iRij, [A,Rij] = iRj, [A,Sij] = 0. (3.34) 

The operators R/ and R are orthogonal to S; in accordance with ()2.48|) . Besides, R/ ■ R = 
3X'^'^'-, and the Casimir operator is /C = Sf + Rf + R^ + = 4. 

Like in the case of double 5*0(4) symmetry studied above, the second step of RG 
procedure generates additional vertices in the exchange part of the interaction Hamilto- 
nian ^i^, 

Hcot = Jl^l ■ Si + J2H1 ■ Si + J3(Ri ■ Sri + 1^2 ' Sir) + ^48^ ■ + J5R ■ 

+ Jei^ll ■ Srl + R2Z ■ Sir) + JjSl ■ {Sir + Sri) + ^8(^-1 ' Sh- + R2 " Sri) 

+ JgRi ■ Sr + JloR ■ S,. + JnRi ■ {Sir + S^) + JuCR-ll " Sir + ^21 ' Srl), (3.35) 

where R^ = X'^''^''Ri, II21 = R2X'^'"^'. The scaling properties of the system are deter- 
mined by a system of 12 scaling equations with initial conditions 

.h{D) = jf, J2{D) = jr, MD)=arJir, MD) = (z = 4 - 12) (3.36) 
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(see Eq. 13. 131 for definitions) specifically 



d\n.d 

dj2 
dliad 

djs 
dlnd 

dji 
dhid 

djb 
dlnd 

dje 
dlnd 
_dh_ 
dhid 

djs 
dlnd 
_dh_ 
d\n.d 

djio 
dlnd 

djn 
dlnd 

dji2 
dlnd 



fi + Ji + jI + Jr + Ju + Jn (J6 + J12) + 
[2(jij2 + Jrjn) + jiUe + ju) - mrj^], 

[jsUi + ji) + jiUb + jio) - mAjQ + ill)], 

J4 + 37 + jI + J?o + Jn + Jn (je + J12) + 
[2jij5 + jji^jz + is) - rnirji], 



jI + jI + jI + Ju 



jI + jI + jI + Ju 



[jeUi + J4) - mirjs], 

(j3 + J8)(j5 + jio) + (j2 + j9)(j6 + ju) 



+ J7(jl + j4) + jll(j2 + jg) 



-b8(ji + J4) + J7(j5 + Jio) - mir{jn + ju)], 
■[2(j4j9 + jrju) + jrUe + ju) - mrjio], 
■[2j4jio + jrUs + is) - mirjg], 
-[jii(ji + J4) + J7(j2 + Jg)], 

■[juUl +j4) - mirjs]- 



(3.37) 



Here the terms proportional to mir arise because the second term in the Hamiltonian 
fl3.30|) contains non-diagonal terms. 

From equations ()3.37|1 . one deduces the Kondo temperature, 

1 

Dil =^ ] V^rnir . (3.38) 



T, 



K2 



Jl + J2 + V(jl+j2)2 + 2j| 



Similarly to the previous case, this equation transforms into the usual exponential form 
when the / and r channels are independent, 

2 

(3.39) 



Tr,. = De ^'1 + -^2 + VO'i + Js)^ + 2j| _ 



lim iK2 

mj^— >0 

Upon increasing mir, the symmetry reduces from SO (5) to 5*0(4). The same happens at 
small mir but with increasing 61 = Es^ —Eti- In the latter case the energy is quenched, 
and at 5/ > Tk2 Eq- (!3.38p transforms into Tk = 6iexp{-[ji{5i) + jsi^i)]^'^} (cf. [73J. 
On the other hand, upon decreasing 6r = E^^ — Esi the symmetry P x 50(4) x 50(4) 
is restored at 5r < Tko- The Kondo effect disappears when 61 changes sign (the ground 
state becomes a singlet). 
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3.2.3 S0{7) Symmetry 



The next asymmetric configuration is illustrated by the flow diagram of Fig lH.fil 
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Figure 3.6: Scaling trajectories for 5*0(7) symmetry in the SW regime. Inset: Zoomed in 
avoided level crossing pattern near the SW line. 

In this case, the manifold {Ti, Si,Tr} is involved in the dynamical symmetry of TQD. 
The relevant symmetry group is 5*0(7). It is generated by six vectors and three scalars. 
These are spin operators Sa (a = /, r) and R-operator R; (see Eq. I3.15|l plus three vector 
operators Rj and three scalar operators Ai involving / — r permutation. Here are the 
expressions for the spherical components of these vectors via Hubbard operators, 

Rf = V2(X^'-°' + X^'~^^), Rl = X^'i"- - X^"-^', 

R^ = V2(X^'°'- + XO'-^O, R2 = - (3-40) 



Ri 



^(^X'^rSi _ x^i^-)^ Rl = -(X^rSl ^ x^l^ 



The scalar operators Ai, A2, A3 now involve the / — r permutations for the triplet states. 
They are defined as 



Ai 

A2 
A, 



2 V 
2 V 



(3.41) 



The (somewhat involved) commutation relations of 07 algebra for these operators and 
various kinematic constraints are presented in Appendix O The SW transformation 
results in the effective Hamiltonian 

A=Ti,Si,Tr ka a=l,r 

+ Jla^a ■ + J2 ^ad " Saa + <^3(R'3 ^ ' ^rl + R 



of ■ Sir) + ■ Si, (3.42) 



a=Lr 



a=Lr 
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where Jia = , J2 = Jir, J3 = otiJir, Ji = otiJf and S„a = X'^'^'^^Sg. The spherical 
components of the vector operators Ri and R2 are 

R'r = iR'sh\ Rr = iR?y, R's = iRSy- (3.43) 

It is easy to see that S^^ + S^; = Ri + R2 and R3 = 'bS^'^ + Ra^"*. 

Like in the case of 5*0(5) symmetry, the tunnehng terms Mj^X^"-^" generate additional 
vertices in the renormalized Hamiltonian Hcot- The number of these vertices and the 
corresponding scahng equations is too wide to be presented here. We leave the description 
of RG procedure for 5*0(7) group for the next chapter (as well as the case of TQD with 
odd occupation), where the case of M^. = is considered. In that situation the scaling 
equations describing the Kondo physics of TQD with 5*0 (n) symmetry are more compact. 

3.3 Conclusions 

The basic physics for all SO{n) symmetries is the same, and we summarize it here. The 
TQD in its ground state cannot be regarded as a simple quantum top in the sense that 
beside its spin operator other vector operators Rn are needed (in order to fully determine 
its quantum states), which have non-zero matrix elements between states of different spin 
multiplets (5*iMj|R„|5*jMj) 7^ 0. These "Runge-Lenz" operators do not appear in the 
isolated dot Hamiltonian (so in some sense they are "hidden"). Yet, they are exposed 
when tunneling between the TQD and leads is switched on. The effective spin Hamiltonian 
which couples the metallic electron spin s with the operators of the TQD then contains 
new exchange terms, J„s ■ R„ beside the ubiquitous ones JjS ■ Sj. The operators Sj and 
R„ generate a dynamical group (usually 5*0 (n)). 

We have analyzed several examples of TQD with even occupation in the parallel 
geometry (Fig. I3.1|l . Our analysis demonstrates the principal features of Kondo effect in 
CQD in comparison with the conventional QD composed of a single well. These examples 
teach us that in Kondo tunneling through CQD, not only the spin rotation but also the 
"Runge-Lenz" type operators R and R are involved. Physically, the operators R describe 
left-right transitions, and different Young schemes give different spin operators in the 
effective co-tunneling Hamiltonians (see Appendix lE]). 
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Chapter 4 
Trimer in Series 



In this Chapter we expose the physics of Kondo tunneling through sandwich-type molecules 
adsorbed on metallic substrate. In Sec. \4.1\ we introduce the system of our study and 
demonstrate the correspondence between the low-energy tunnel spectra of chemis orbed lan- 
thanocene molecule and an artificial trimer, i.e., TQD in series geometry. In Sec. \4-^ 
we concentrate on the case of even occupation. The scaling equations are derived and 
the Kondo temperatures are calculated for the evenly occupied trimer in the cases of the 
P X ^0(4) X SO{A), 50(5), soil) and P x 50(3) x 50(3) dynamical symmetries. The 
dynamical- symmetry phase diagram is displayed and the possibility of its experimental re- 
alization is outlined. The anisotropic Kondo effect induced by an external magnetic field 
is discussed in Sec. \4-3[ It is shown that the symmetry group for such magnetic field 
induced Kondo tunneling is SU{3). The case of odd occupation is considered in Sec. \4-4\ 
When the ground state of the trimer is a doublet, the effective spin Hamiltonian of the 
trimer manifests a two-channel Kondo problem albeit only in the weak coupling regime. 
Analysis of the Kondo effect in cases of higher spin degeneracy of the trimer ground state 
is carried out in relation with dynamical symmetries. In the Conclusions we underscore 
the main results obtained. 

4.1 Introduction 

In this Chapter we extend the theory of single-electron tunneling developed in Chapter El 
for complex quantum dots to the case of sandwich-type molecules adsorbed on metallic 
substrate. The geometry of the nano-objects under investigation is: metallic subsrate 
(MS) - molecule - nanotip of scanning tunnel microscope (STM). The specific objects of 
our studies are lanthanocene molecules Ln(C8H8)2 where the magnetic ion Ln=Ce, Yb 
(central ion) is secluded in a cage of carbon-containing radicals and only these radicals are 
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in direct tunnel contact with MS and STM. Ce is known as a mixed-valent ion in many 
molecular and crystalline configurations. This means that the covalent bonding, i.e., 
hybridization between strongly correlated 4f-electron and weakly interacting molecular 
vr-orbitals Mo = CgHg in cerocene molecule Ce(C8H8)2 characterized by the parameter 
Vh, is noticeable. Besides, the difference between the ionization energy £/ of 4f-electron 
and the ionization energy of an electron in molecular vr-orbital is large enough, so that 
V/i <^ — £/. As is shown in jlSl IIHI 1^0], the direct consequence of this inequality is 
unusual electronic and spin structure of the ground state and low-energy excitations in 
Ce(C8H8)2- Instead of purely ionic bonding Ce^+(4f^)(Mo)2~(e2„) a mixed valence state 
arises with admixture of configuration 4f''e2„ to the ground state spin singlet. Due to the 
Pauli principle, such admixture is forbidden in a triplet state. As a result, the energy 
difference between the ground state singlet and excited triplet is controlled by the small 
parameter Vh/{e-„ — e j). 

This scenario reminds the mechanism of singlet-triplet splitting in asymmetric DQD 
(Sec. 12. 3|) and exactly coincides with that for TQD with small central dot, whose prop- 
erties are studied in Chapter 3. The singly occupied central dot with strong Coulomb 
blockade plays the same part as magnetic ion (Ce or Yb) in lantanocene molecules. Two 
side dots in TQD play the same role as two molecular rings CsHg in formation of low 
energy spin spectrum. Only a pair of the highest occupied molecular orbitals (HOMO) 
is involved in formation of this spectrum. All other states may be treated as molecular 
excitons, which are quenched within the scale of Vj^/^e.,, — Sf) above the ground state. As 
a result, the orbital degrees of freedom are irrelevant in the Kondo regime. The similarity 
between the configurations "STM-Ce(C8H8)2-metaUic adsorbent" and source - TQD - 
drain" is illustrated by Figs. I4.1l and l4.2[ The only essential condition for modelling both 
systems by the same Hamiltonian ()4.1|) and treating them as a trimer in contact with 
metallic leads is the inequality W > V, i.e., the demand that the tunnel contact with the 
reservoir does not destroy the coherent quantum mechanical state of a trimer. 

It was mentioned already in Chapter IHl that a TQD with leads / and r representing 
independent tunneling channels can be mapped onto a TQD in a series by means of 
geometrical conformal transformation. Indeed, if the inter-channel tunneling amplitude 
tir in the Hamiltonian ()3.3j) vanishes, one may apply a rotation in source-drain space 
separately to each channel and exclude the odd s — d combination of lead states both in 
the /- and r-channel j^E]- Since now each lead is coupled to its own reservoir, and one 
arrives at the series configuration shown in Fig. 14.21 

It is virtually impossible to conceive an additional transformation after which the odd 
combination of lead states is excluded from the tunneling Hamiltonian [75^. As a result. 
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Figure 4.1: A molecule with strong correlations is modelled by a TQD in a series geometry. 

the challenging situation arises in case of odd occupation A/" = 3, where the net spin 
of TQD is S = 1/2, and the two leads play part of two channels in Kondo tunneling 
Hamiltonian. Unfortunately, despite the occurrence of two electron channels in the spin 
Hamiltonian, the complete mapping on the two-channel Kondo problem is not attained 
because there is an additional cotunneling term JirS-sir + H.c. (s;^ is determined by (jSHEI)) 
which turns out to be relevant, and the two-channel fixed point cannot be reached (see 
Sec. 14. 4j) . And yet, from the point of view of dynamical symmetry the series geometry 
offers a new perspective which we analyze in the present chapter for the cases of even and 
odd occupation. 



4.2 Even Occupation 

Consider then a trimer in series (Fig. 14. 2j) with four electron occupation A/" = 4. The 
Hamiltonian of the system can be written in the form, 

Aa A k<T b=s,d 

+ T.T.i(^"^''^tk. + Vr^''c^,jX'^ + H.c.]. (4.1) 

AA kcr 

Here |A), |A) are the four- and three-electron eigenfunctions ()A.5|1 and ()A.9|1 . respectively; 
i?A) Ex are the four- and three-electron energy levels, respectively (five electron states 
cost much energy and are discarded); X^^ = |A)(A| are number changing dot Hubbard 
operators. The tunneling amplitudes V^J^ = Va{X\daa\-^) {a = l,r) depend exphcitly on 
the respective 3 — 4 particle quantum numbers A, A. Note that direct tunneling through 
the TQD is suppressed due to electron level mismatch and Coulomb blockade, so that 
only cotunneling mechanism contributes to the current. 
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Figure 4.2: Triple quantum dot in series. Left (/) and right (r) dots are coupled by 
tunneling Wi^r to the central (c) dot and by tunneling Vi^r to the source (s) (left) and 
drain (d) (right) leads. 

After a SW transformation the generic Hamiltonian ()3.12p simplifies in this case to 

Aa kcr b=l,r a=l,r 

+ ^ Ja^^a ■ Sa + Jlr ^ R-a " Saa, (4.2) 
a=l,r a=l,r a=l,r 

(the notation l,r is used for the electron states both in the leads and in the TQD). The 
antiferromagnetic coupling constants are defined by (jSHSl)- The vectors Sa, Rq and Ra 
are the dot operators ()3.15p . P is the permutation operator ()3.14|) . and the components 
of the vectors s^, Saa are determined in Eqs. ()2.(i3j) (with a = a = l,r) and ()3.1fi|l . The 
vector operators S^, R^, Rq and the permutation operator P manifest the dynamical 
symmetry of the TQD. 

We now discuss possible realization of P x 5*0(4) x 50(4), 50(5), 50(7) and P x 
5*0(3) X 5*0(3) symmetries arising in the TQD with A/" = 4 [ZH El . Due to the absence of 
interchannel mixing, the avoided crossing effect does not arise in the series geometry. As 
a result, the cases of P x 50(4) x 50(4), 50(5) and 50(7) symmetry are characterized 
by the same flow diagram of Figs. 13. 3| 13.51 and 13. (il but without avoided crossing effects 
shown in the insets. 

Let us commence the analysis of the Kondo effect in the series geometry with the case 
P X 5*0(4) X 5*0(4) where Ei^ = Ere and Ft, = Fr, (Fig. EISl)- In this case the exchange 
part of the Hamiltonian ()4.2|1 is a simplified version of the Hamiltonian ()3.23j) with the 
boundary conditions ()3.24|) . The scaling equations are the same as ()3.25|) with rriir = 0. 
Solving them one gets Eq. ()3.27j) for the Kondo temperature. 
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When ~ -Et; ~ < E^^ (Fig. 13. 5|) . the TQD possesses the S0{5) symmetry. In 
this case the interaction Hamiltonian has the form 



Hcot = Jl^l ■ Si + J2Rz ■ S; + J3(Ri ■ + R2 ■ Sir 



(4.3) 



which is the same as in Eq. ()3.30p with Ri, R2 determined by Eq. ()3.3H) . Respectively, the 
effective Hamiltonian for the Anderson scaling is a reduced version of the Hamiltonian 



H, 



cot 



JlSi ■ Si + J2R1 ■ Si + J3(Ri ■ Sri + R2 ■ Sir) + -^48^ ■ S, 



(4.4) 



with the boundary conditions ()3.36|1 for Jj, i = 1 — A. 
The scaling equations have the form 



dji 



dlnd 

dj2 
dlnd 

dj3 
dlnd 

dhid 



-2jlj2, 
-J3(jl + j4), 



(4,5) 



Of course, Eqs. ()4.5j) for the Kondo temperature yield the limiting value ()3.39|1 . 

When ^ Et, ^ Es, < Es, (Fig. EIHl), the TQD possesses the SO (7) symmetry. In 
this case the Anderson RG procedure adds three additional vertices in the exchange part 
of the basic SW Hamiltonian ()3.42|) . 



Hcot — Jla^a ■ Sa + </2 ^aa ' ^aa + <^3(R3 ^ ' + R3 '' ' Sir 

a=l,r a=l,r 

+ JiRl ■ Si + JsaS + JeRi ■ Sr. 



(4.6) 



a=Lr 



The boundary conditions for solving the scaling equations are 

Jla{D) = jJ, MD) = Jir, J3{D) = aiJlr, 



(4.7) 
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The system of scaling equations 
djii 



dlnd 
djir 
dlnd 

_dh_ 
dlnd 

djs 
dlnd 

djA 
dlnd 

djbi 
dlnd 

djbr 

dlnd 

dje 
dlnd 



j2Ull + jlr + jbl + jbr) + hiii + 3%) 



[j2(j4 + Je) + JsO'lr + J5r)] , 



- [2jl/j4 + 32h] 



ill 



■2 

J5r + 2 + 2 



[2j5/j6 + J2J3] 



(4. 



is now solvable analytically, and the Kondo temperature is, 

4 



K 



D exp 



(4.9) 



where j+ = ju + J4 + jir, j_ = ji; + J4 - jir- 

Like in the cases considered above, the Kondo temperature and the dynamical sym- 
metry itself depend on the level splitting. On quenching the Si state (increasing 6ir = 

— Et^), the pattern is changed into a P x 5*0(3) x 5*0(3) symmetry of two degenerate 
triplets with a mirror reflection axis. Changing the sign of 6ir one arrives at a singlet 
regime with = 0. 

When the lowest renormalized states in the SW limit are two triplets Ti and T^, 
the TQD possesses the P x 5*0(3) x 5*0(3) symmetry with mirror reflection axis. The 
corresponding co-tunneling spin Hamiltonian has the form. 



a=l,r 



(4.10) 



a=Lr 



Here J^a{D) = KV(eF - Ca), ^2(5) = ^Eai^F - ea)'' and Jsa{D) = 0. The P x 
5*0(3) X 50(3) symmetry is generated by the spin one operators Sa with projections 
fia = la, Oq, Iq, and the left-right permutation operator P ()3.14|) . 
The system of scaling equations for the Hamiltonian ()4.10|1 is. 
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Figure 4.3: Phase diagram of TQD. The numerous dynamical symmetries of a TQD in 
the parallel geometry are presented in the plane of experimentally tunable parameters 
X = Ti/Tr and y = Eid Ere- 

where j = poJ, a = l,r. From Eqs. ()4.1ip we obtain the Kondo temperature, provided 
\Eti — Et^\ < Tk, 



The results of calculations described in this section are summarized in Fig. 14.31 The 
central domain of size Tko describes the fully symmetric state where there is left-right 
symmetry. Other regimes of Kondo tunneling correspond to lines or segments in the {x, y} 
plane. These lines correspond to cases of higher conductance (ZBA). On the other hand, at 
some hatched regions, the TQD has a singlet ground state and the Kondo effect is absent. 
These are marked by the vertically hatched domain. Both the tunneling rates which 
enter the ratio x and the relative level positions which determine the parameter y depend 
on the applied potentials, so the phase diagram presented in Fig. 14.31 can be scanned 
experimentally by appropriate variations of and Vga- This is a rare occasion where an 
abstract concept like dynamical symmetry can be felt and tuned by experimentalists. The 
quantity that is measured in tunneling experiments is the zero-bias anomaly (ZBA) in 
tunnel conductance g [23 123- The ZBA peak is strongly temperature dependent, and 
this dependence is scaled by Tk- In particular, in a high temperature region T > Tk, 
where the scaling approach is valid, the conductance behaves as 



As it has been demonstrated above, Tk in CQD is a non-universal quantity due to partial 
break-down of dynamical symmetry in these quantum dots. It has a maximum value in 
the point of highest symmetry P x SO (4) x SO {4), and depends on the parameters 6a 



Tk = D exp 



2 



(4.12) 



3ii + jiT + - jivY + 2j| 



g{T)^\Ti-\T/TK). 



(4.13) 
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in the less symmetric phases (see, e.g., Eqs. 13.271 13.391 14. 12^ . Thus, scanning 
the phase diagram means changing TxiSa)- These changes are shown in Fig. 14.41 which 



Tk/Tko 




Figure 4.4: Variation of Kondo temperature with 6ri = Vgr — Vgi. Increasing this param- 
eter removes some of the degeneracy and either "breaks" or reduces the corresponding 
dynamical symmetry. 

illustrates the evolution of Tk with 6ri for x = 0.96, 0.8 and 0.7 corresponding to a 
symmetry change from P x 50(4) x 50(4), 50(7) to P x 50(3) x 50(3) and from 
50(5) to 50(4), respectively. It is clear that the conductance measured at given T 
should follow variation of Tk in accordance with ()4.13j) . 

4.3 Anisotropic Kondo Tunneling through Trimer 
4.3.1 Generalities 

In all examples of CQDs considered above the co-tunneling problem is mapped on the 
specific spin Hamiltonian where both S and R vectors are involved in resonance cotun- 
neling. There are, however, more exotic situations where the effective spin Hamiltonian is 
in fact a "Runge-Lenz" Hamiltonian in the sense that the vectors R alone are responsible 
for Kondo effect. Actually, just this aspect of dynamical symmetry in Kondo tunneling 
was considered in the theoretical papers ^21 ^1 ^1 QB] and observed experimentally in 
Refs. [T^ 1^2] , in which the Kondo effect in planar and vertical QDs induced by external 
magnetic field B has been studied. In this section we lay down the theoretical basis for 
this somewhat unusual kind of Kondo effect. 

Consider again the case of TQD in series geometry with A/" = 4. In the previous 
sections the variation of spin symmetry was due to the interplay of two contributions to 
indirect exchange coupling between the spins Sa- One source of such an exchange is tun- 
neling within the CQD (amplitudes Wa) and another one is the tunneling between the dots 
and the leads (amplitudes Va). An appropriate tuning of these two contributions results 
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in accidental degeneracy of spin states (elimination of exchange splitting), and various 
combinations of these accidental degeneracies lead to the rich phase diagram presented in 
Fig. 14.31 A somewhat more crude approach, yet more compatible with experimental ob- 
servation of such interplay is provided by the Zeeman effect. This mechanism is effective 
for CQD which remains in a singlet ground state after all exchange renormalizations have 
taken place. The negative exchange energy 6a may then be compensated by the Zeeman 
splitting of the nearest triplet states, and Kondo effect arises once this compensation is 
complete jT^j. From the point of view of dynamical symmetry, the degeneracy induced 
by magnetic field means realization of one possible subgroup of the non-compact group 
SO{n) (see Eq. 12.641 and corresponding discussion in Section \2A\\ . The transformation 
5*0(4) — > SU{2) for DQD in magnetic field was discussed in Ref. [SE]- 

4.3.2 Trimer with SU(3) Dynamical Symmetry 

In similarity with DQD, the Kondo tunneling may be induced by external field B in the 
non-magnetic sector of the phase diagram of Fig. 14.31 A very peculiar Kondo tunneling is 
induced by an external magnetic field B in the non-magnetic sector of the phase diagram 
of Fig. 14.31 close to the 50(5) line. In this case, a remarkable symmetry reduction occurs 
when the Zeeman splitting compensates negative 5i^r = ^ — ■ Then we are left with 
the subspace of states {Tl;, S";, S'r}, and the interaction Hamiltonian has the form, 

+ Ji{R^sl + Risli) + {J^Rl + J(,Rl)sl + J7{Rts; + R{st). (4.14) 

Here 

MD) = MD) = ^, MD) = jr, 

MD) = arJir, MD) = (« = 4-7). (4.15) 

The operators Ri, R2, -R3 and R4 are defined as, 

Rl = ^{X^'^' - X'^'^'), i?+ = Xi'^S R^ = {R+y, 
Rl = ^(X^'i' - X^^^O, Rt = X'-'^, R^ = (i?+)t, 

i?3 = i?4 = ^x^'-^'. (4.16) 

We see that the anisotropic Kondo Hamiltonian ()4.14|) is quite unconventional. There are 
several different terms responsible for transverse and longitudinal exchange involving the 
R-operators which generate both S^/T and S^/Sa transitions. 
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The operators ()4.16|) obey the following commutation relations, 



[Rlj,Rlk] — iejkrnRlm, [R2j,R2k\ — i^jkmR2mi 



[Rij, R2k] = -^{Rs — RijSjkil — Sjz 



+ -Gjkm(^Rlm^kz + R2m^jz —^mz{R-i + -^4)^ ; 

1 ^ c- 



[Rljj -R3] = —T^R^Sjz + —riR2x + 'i'R2y){.5jx — iSjy), 



[Rlj, -R4] = -^Ri^jz —{R2X — iR2y)iSjx + iSjy), 

1 ^/3 
[i?2i, R3] = ^^S^jz -^{Rlx + iRly){Sjx + i6jy), 

1 

[i?2j, R4] = "^Ri^jz + —^{Rlx + iRly){Sjx — i^jy), 

[Rs,R,] = l{R',-Rl). (4.17) 

These operators generate the algebra U3 in the reduced spin space {Tli, Si, Sr} specified 
by the Casimir operator 

3 

+ R2 + -R3 + -R4 = - . 

Therefore, in this case the TQD possesses SU (3) symmetry. These R operators may be 



represented via the familiar Gell-Mann matrices Aj 


(z = l,.. 


.,8) for the SU{3) group 


Rt = ^(Ai + ^s), 


Ri 




- ^Aa) , 


vz A3 

^1 - y, 


i?2 = 


= + 


VSAg), 


Rt = ^(A4 + iA5), 


R2 




- ^As) , 


R3 = ^{K + ^X7), 


R4 = 




! - ^At) . 



As far as the RG procedure for the "Runge-Lenz" exchange Hamiltonian ()4.14p is 
concerned, the poor-man scaling procedure is applicable also for the R operators. The 
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scaling equations have the form 
dji 



d\nd 
dhid 

d\n.d 
dji 

d\n.d 
djb 

dlnd 
djr 

dlnd 



dj2 



dlnd 



J3l + ^) ^ —JSrJi 



J3r [-y^Sl + ^J7 



(ilnc? 



J5J7 + 7rJ3rj4 



(4.18) 



where j = PqJ. We cannot demonstrate analytical solution of this system, but the 
numerical solution shows that stable infinite fixed point exists in this case like in all 
previous configurations. 

Another type of field induced Kondo effect is realized in the symmetric case of 5 = 
— „ < 0. Now the Zeeman splitting compensates negative 6. Then the two 
components of the triplets, namely -^Tig „ cross with the singlet state energy , and the 
symmetry group of the TQD in magnetic field is SU{3) as in the case considered above. 



4.3.3 Summary 

It has been demonstrated that the loss of rotational invariance in external magnetic field 
radically changes the dynamical symmetry of TQD. We considered here two examples of 
symmetry reduction, namely S0{5) SU{3) and P x ^0(4) x ^0(4) -> SU{3). In all 
cases the Kondo exchange is anisotropic, which, of course, reflects the axial anisotropy 
induced by the external field. These examples as well as the SO (4) SU{2) reduction 
considered earlier [321 EH] describe the magnetic field induced Kondo effect owing to the 
dynamical symmetry of complex quantum dots. Similar reduction SO{n) SU{n') 
induced by magnetic field may arise also in more complicated configurations, and in 
particular in the parallel geometry. The immense complexity of scaling procedure adds 
nothing new to the general pattern of the field induced anisotropy of Kondo tunneling, 
so we confine ourselves with these two examples. 

Although the anisotropic Kondo Hamiltonian was introduced formally at the early 
stage of Kondo physics [1011^, it was rather difficult to perceive how such Hamiltonian 
is derivable from the generic Anderson-type Hamiltonian. It was found that the effec- 
tive anisotropy arises in cases where the pseudo-spin degrees of freedom (like a two-level 
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system) are responsible for anomalous scattering. Another possibility is the introduc- 
tion of magnetic anisotropy in the generic spin Hamiltonian due to spin-orbit interaction 
(see Ref. jH2] for a review of such models). One should also mention the remarkable 
possibility of magnetic field induced anisotropic Kondo effect on a magnetic impurity in 
ferromagnetic rare-earth metals with easy plane magnetic anisotropy This model is 
close to our model from the point of view of effective spin Hamiltonian, but the sources 
of anisotropy are different in the two systems. In our case the interplay between singlet 
and triplet components of spin multiplet is an eventual source both of the Kondo effect 
itself and of its anisotropy in external magnetic field. Previously, the manifestation of 
SU{3) symmetry in anisotropic magnetic systems were established in Refs. jSH I85j. It 
was shown, in particular, that this dynamical symmetry predetermines the properties of 
collective excitations in anisotropic Heisenberg ferromagnet. In the presence of single-ion 
anisotropy the relation between the Hubbard operators for S* = 1 and Gell-Mann matrices 
A were established. It worth also mentioning in this context the SU{A) D 50(5) algebraic 
structure of superconducting and antiferromagnetic coherent states in cuprate High-Tc 
materials |86j . 



4.4 Odd Occupation 

We now turn our attention to investigation of the dynamical symmetries of TQD in series 

with odd occupation Af = 3, whose low-energy spin multiplet contains two spin 1/2 

doublets I -81,2) and a spin quartet \Q) with corresponding energies 

3 

Eb, = Ec + ei+Er- -[Wl(3l + Wr(3r], 

Eb, = 6, + 6i+6r-^[Wif3i + Wr(3r], 
Eq = Ec + ei+Er. (4.19) 

There are also four charge-transfer excitonic counterparts of the spin doublets separated 
by the charge transfer gaps ei — Ec + Qi and Sr — Ec + Qr from the above states (see 
Appendix 

Like in the four-electron case, the scaling equations ()2.33|1 may be derived with different 
tunneling rates for different spin states {Tq for the quartet and T^. (i = 1,2) for the 
doublets). 

Tq = npo {v^' + K') , r^, = 7?rQ, r^, = ^ITq, (4.20) 

with 

7i = yi^|c^ + ^, I2 = ^l-l{f3f + m- (4.21) 
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Since Tq > F^j, F^j, the scaling trajectories cross in a unique manner: This is the com- 
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Figure 4.5: Scahng trajectories resulting in 50(4) x 5?/ (2) symmetry of TQD with A/" = 3. 

plete degenerate configuration where all three phase trajectories intersect \Eq{D*^ = 
Ebi{D*) = Eb2{D*)] at the same point D*. This happens at bandwidth D = D* iVw WJ^ 
whose value is estimated as 



D* = Doexp 



vrr 



(4.22) 



where 



W^Erc+W^Ei, 



ElcErc- 



W^El + W^El 

This level crossing may occur either before or after reaching the SW limit D where scaling 
terminates Below we discuss the Kondo physics arising in the cases: D* < D, D* = D 
and D* > D. 



4.4.1 Towards Two-Channel Kondo Effect 

When D* < D, the lowest renormalized state in the SW limit is a doublet Bi. Following 
an RG procedure and a SW transformation, the spin Hamiltonian in this case reads 

Hcot = <^/S ■ S/ + JrS ■ Sr + Jlr-S ■ (s^^ + Sri). (4.23) 

The spin 1/2 operator S acts on \Bi„=-^^i) fEq. ()A.9|) ). whereas the lead electrons spin 
operators and Saa are determined in Eqs. ()2.fi3j) (with a = a = l,r) and ()3.1(')|) . The 
exchange coupling constants are 

J.^-^, J^.^-f^. (4.24) 

3(eF - So) ^{ep - So) 

The Hamiltonian ()4.23j) then encodes a two-channel Kondo physics, where the leads serve 
as two independent channels and Tr- = max{TKu Trt} with T^a = -De"^/-'" and = poJa- 
A poor-man scaling technique is used to renormalize the exchange constants by reduc- 
ing the band- width D ^ D. The pertinent fixed points are then identified as D Tk 
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Unlike the situation encountered in the single- channel Kondo effect, third order dia- 
grams in addition to the usual single-loop ones should be included (see Fig. 5 in Ref. [66j 
and Fig. 9 in Ref. [HZl)- With a = l,r and a = r,l the three RG equations for ji,jr,jir 
are 



dja 

dlnd 
djir 
dlnd 



-Jlr (jl + Jr) + JlrUf + Jr + '^jfr) ■ 



(4.25) 



On the symmetry plane ji = jr = j, Eqs. ()4.25|) reduce to a couple of RG equations for 

jl,2 = j ± jlr 

dji 



d\n.d 



1,2), 



(4.26) 



subject to ji{D = D) = jiQ = poJi. These are the well-known equations for the anisotropic 
two-channel Kondo effect [HE]- With 0j = (ji + j 2 — l)/ji) Ci = 0jo ^ 4>io {(pio = (pi{D)) 
and Li{x) = x — ln(l + Ci/x) — 21na;, the solution of the system ()4.26|) is 



Li((f)j) - Li 



ho 



) = In 



1,2). 



(4.27) 



The scaling trajectories in the sector {ji > j^ > 0,jir = 0) and in the symmetry plane 
with < jlr < j are shown in Fig. 14.61 Although the fixed point (1/2, 1/2, 0) remains 




0.5 1 1.5 2 

Figure 4.6: Scaling trajectories for two-channel Kondo effect in TQD. 

inaccessible if jir 7^ 0, one may approach it close enough starting from an initial condition 
jiro jio,jro- Realization of this inequality is a generic property of TQD in series shown 
in Fig. 14.21 A similar scenario was offered in [HHl IHO] for a QD between two interacting 
wires. 

According to general perturbative expression for the dot conductance [HH] , its zero-bias 
anomaly is encoded in the third order term. 



G^'^ = Gojfr[jiiT)+jr{T)], (Go = ^) 



(4.28) 
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Here the temperature T replaces the bandwidth D in the solution ()4.27|) . Let us present 
a qualitative discussion of the conductance G[ia{T)] (or in an experimentalist friendly 
form, G{vga,T)) based on the flow diagram displayed in Fig. 14.61 (Strictly speaking, 
the RG method and hence the discussion below, is mostly reliable in the weak-coupling 
regime T > Tk)- Varying T implies moving on a curve [ji{T), jr{T), jir{T)] in three- 
dimensional parameter space (Fig. 14. 6|) . and the corresponding values of the exchange 
parameters determine the conductance according to equation ()4.28j) . Note that if, initially, 
jio = jro = jo the point will remain on a curve [j(T), j(T), jir{T)] located on the symmetry 
plane. By varying Vga it is possible to tune the initial condition {jio,jro) from the highly 
asymmetric case jio ^ jro to the fully symmetric case jio = jro- For a fixed value of jiro 
the conductance shoots up (logarithmically) at a certain temperature T* which decreases 
toward with |j;o ~ jro\ and j/^o- The closer is T* to T^, the closer is the behavior 
of the conductance to that expected in a generic two-channel situation. Thus, although 
the isotropic two-channel Kondo physics is unachievable in the strong-coupling limit, its 
precursor might show up in the intermediate-coupling regime. 




5 10 15 20 



Figure 4.7: Conductance G in units of Go as a function of temperature (r = T/Tk), at 
various gate voltages. The lines correspond to: (a) the symmetric case ji = jr {vgi = Vgr), 
(b-d) ji 3> jr, with Vgi — Vgr = 0.03, 0.06 and 0.09. At r — cx) all lines converge to the 
bare conductance. 

The conductance G{vga,T) as function of T for several values of Vga and the same 
value of jiro is displayed by the family of curves in Fig. 14.71 For G displayed in curve 
a, T* /Tk ~ 3 and for T > T* it is very similar to what is expected in an isotropic 
two-channel system. Alternatively, holding T and changing gate voltages Vga enables 
an experimentalist to virtually cross the symmetry plane. This is equivalent to moving 
vertically downward on Fig. 14.71 At high temperature the curves almost coalesce and 
the conductance is virtually fiat. At low temperature (still above Tk) the conductance 
exhibits a sharp minimum. This is summarized in Fig. 14.81 
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Figure 4.8: Conductance G in units of Go as a function of gate voltage at various tem- 
peratures (at the origin ji = jr). 

4.4.2 Higher Degeneracy and Dynamical Symmetries 

If the degenerate point ()4.22j) occurs in the SW crossover region, i.e., if D* ^ D, the 
SW procedure involves all three spin states (Fig. I4.5p . and it results in the following 
cotunneling Hamiltonian 

Hcot = Yl (-Ja^ + -Ja^^) ■ (4-29) 
a=l,r 

where S is the spin 1 operator and R is the R-operator describing S/T transition similar 
to that for spin rotator [36 . The coupling constants are 

Jl = ^7r^y Ja^ = ^^Jlr- (4.30) 

This is a somewhat unexpected situation where Kondo tunneling in a quantum dot 
with odd occupation demonstrates the exchange Hamiltonian of a quantum dot with even 
occupation. The reason for this scenario is the specific structure of the wave function 
of TQD with A/" = 3. The corresponding wave functions |A) (see Appendix ^ are 
vector sums of states composed of a "passive" electron sitting in the central dot and 
singlet /triplet (S/T) two-electron states in the l,r dots. Constructing the eigenstates |A) 
using certain Young tableaux (see Appendix D), one concludes that the spin dynamics 
of such TQD is represented by the spin 1 operator S corresponding to the / — r triplet, 
the corresponding R-operator R and the spin 1/2 operator Sc of a passive electron in the 
central well. The latter does not enter the effective Hamiltonian Hcot ()4.29|) but influences 
the kinematic constraint via Casimir operator /C = + + = ^. The dynamical 
symmetry is therefore SO (4) x SU{2), and only the SO (4) subgroup is involved in Kondo 
tunneling. 



The scaling equations have the form, 

djla 

dlnd 

dj2a 

dlnd 



b'la + i2a] 5 

-2jlaj2a, (4.31) 
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where jia = PqJI, j2a = Po^a^ (ct = l,r). From Eqs. ()4.3ip we obtain the Kondo 
temperature, 

Tk = max{TKi, Tkt}, (4.32) 

with Txa = I) exp [-1/ + jsa)]- 

An additional dynamical symmetry arises in the case when D* > D. In this case 
the ground state of TQD is a quartet S=3/2, and we arrive at a standard underscreened 
Kondo effect for SU{2) quantum dot as an ultimate limit of the above highly degenerate 
state. 

4.4.3 Summary 

To conclude this section, it might be useful here to underscore the following points: (1) In 
a TQD (Fig. 14. 2|) . the two-channel (left-right leads) Kondo Hamiltonian ()4.23|) emerges 
in which the impurity is a real spin and the current is due solely to co-tunneling. The 
corresponding exchange constant Jir is a relevant parameter: by taking even and odd 
combinations, the system is mapped on an anisotropic two-channel Kondo problem where 
Jir determines the degree of anisotropy. (2) Although the generic two-channel Kondo 
fixed-point is not achievable in the strong coupling limit, inspecting the conductance 
G{vga,T) as function of temperature (Fig. 14. 7|) and gate voltage (Fig. 14. 8|) suggests 
an experimentally controllable detection of its precursor in the weak and intermediate 
coupling regimes. Apparently, genuine multichannel Kondo regime with finite fixed point 
may be achieved for configurations with more than two terminals. (3) There exists a 
scenario of level degeneracy in which TQD with half-integer spin behaves as a dot with 
integer spin in Kondo tunneling regime. 

4.5 Conclusions 

We have analyzed the occurrence of dynamical symmetries in collective phenomena, which 
accompany quantum tunneling through chemisorbed sandwich-type molecules and com- 
plex quantum dots in configurations having a form of linear trimer. These symmetries 
emerge when the trimer is coupled with metallic electrodes under the conditions of strong 
Coulomb blockade in one of its three constituents and nearly degenerate low energy spin 
spectrum. 

As a prototype of complex molecules, where the magnetic ion is sandwiched between 
two molecular radicals, we have chosen lanthanocene molecule Ln(C8H8)2 (Ln=Ce, Yb), 
following P. Fulde's proposal jHO]. This molecule is characterized by anomalously soft 
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spectrum of singlet-triplet excitations. Artificial "Fulde molecule" in a form of linear 
triple dot (Fig. I4.2|l with even electron occupation mimics the low-energy spectrum of 
cerocene. Fulde et. al. |1HI HHI considered the interplay between two spin states (singlet 
and triplet) and singlet charge exciton as a predecessor of genuine Kondo singlet /triplet 
pairs, which arise in classical Kondo effect for a local spin immersed into a Fermi sea. We 
have shown that the Fulde trimer as a whole may be attached to a Fermi basin (metallic 
layer or electrodes). As a result, the pseudo Kondo S/T pairs becomes a new source 
of Kondo scattering/tunneling with quite sophisticated structure of a "scatterer". This 
structure is elegantly described in terms of non-compact dynamical symmetry groups with 
complicated o„ algebras. Application of an external magnetic field results in additional 
accidental degeneracies. Due to the loss of spin-rotational symmetry the effective cotun- 
neling Hamiltonian acquires spin anisotropy, and the dynamical symmetry of a trimer 
is radically changed. Although the main focus in this Chapter is related to the study 
of triple quantum dots, the generalization to other quantum dot structures is indeed 
straightforward. 

The difference between series and parallel geometries of TQD coupled to the leads by 
two channels exists only at non-zero interchannel mixing in the leads, tir 7^ 0. One may 
control the dynamical symmetry of Kondo tunneling through TQD by varying the gate 
voltage and/or lead-dot tunneling rate. In the case of odd electron occupation [M = 3) 
when the ground-state of the isolated TQD is a doublet and higher spin excitations can 
be neglected, the effective low-energy Hamiltonian of a TQD in series manifests a two- 
channel Kondo problem albeit only in the weak coupling regime [ZSj- To describe the 
flow diagram in this case, one should go beyond the one-loop approximation in RG flow 
equations [021 • The nominal spin of CQD does not necessarily coincide with that involved 
in Kondo tunneling. A simple albeit striking realization of this scenario in this context is 
the case of TQD with A/" = 3, which manifests itself as a dot with integer or half- integer 
spin (depending on gate voltages). 

Since we were interested in the symmetry aspect of Kondo tunneling Hamiltonian, we 
restricted ourselves by derivation of RG flow equations and solving them for obtaining 
the Kondo temperature. In all cases the TQDs possess strong coupling fixed point char- 
acteristic for spin 1/2 and/or spin 1 case. We did not calculate the tunnel conductance 
in details, because it reproduces the main features of Kondo-type zero bias anomalies 
studied extensively by many authors (see, e.g., O CH CHI CEl EHl ED] ) • The novel feature 
is the possibility of changing Tx by scanning the phase diagram of Fig. 14.31 Then the 
zero bias anomaly follows all symmetry crossovers induced by experimentally tunable gate 
voltages and tunneling rates. 
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The main message of this Chapter is that symmetry enters the realm of mesoscopic 
physics in a rather non-trivial manner. Dynamical symmetry in this context is not just a 
geometrical concept but, rather, intimately related with the physics of strong correlations 
and exchange interactions. The relation with other branches of physics makes it even 

more attractive. The groups SO{n) play an important role in Particle Physics as well as 
in model building for high temperature superconductivity (especially 5*0(5)). The role 
of the group SU{3) in Particle Physics cannot be overestimated and its role in Nuclear 
Physics in relation with the interacting Boson model is well recognized. This work extends 
the role of these Lie groups in Condensed Matter Physics. 
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Chapter 5 

Kondo Tunneling through Triangular 
Trimer in Ring Geometry 

In this Chapter we consider a ring-like triangular trimer, i.e., triangular triple quantum 
dot (TTQD), and focus on its symmetry properties, which influence the Kondo tunneling. 
The basic concepts are introduced in Sec. I5.il In Sec. \5.'J we construct the Hamiltonian 
of the TTQD both in three- and two-terminal geometry and expose its energy spectrum. In 
Sec. \5.'A we focus on the modification of the symmetry of TTQD in an external magnetic 
field. It is shown that TTQD in a magnetic field demonstrates unique combination of 
Kondo and Aharonov-Bohm features. The poor-man scaling equations are solved and the 
Kondo temperatures are calculated for the cases of SU{2) and SU{4) symmetries. The 
conductance of a TTQD strongly depends on the underlying dynamical symmetry group. 
We show that the interplay between continuous spin-rotational symmetry SU{2), gauge 
symmetry U{1) and discrete symmetry C^^ of triangle may result in sharp enhancement 
or complete suppression of tunnel conductance as a function of magnetic flux through the 
TTQD. In Conclusions we summarize the results obtained. 

5.1 Introduction 

In the previous chapters we studied tunnehng through hnear trimers in parallel and serial 
configurations, where the dots are ordered linearly either parallel or perpendicular to 
the metallic leads. Meanwhile, modern experimental methods allow also fabrication of 
quantum dots in a ring geometry. This ring may have a form of closed gutter pT| W]\ 
or be composed by several separate dots coupled by tunnel channels. In the latter case 
the simplest configuration is a triangle. Triangular triple quantum dot (TTQD) was 
considered theoretically [HSl and realized experimentally very recently |94,, in order to 
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demonstrate the ratchet effect in single electron tunneling. To achieve this effect the 
authors proposed a configuration, where two of the three puddles are coupled in series 
with the leads (source and drain), while the third one has a tunnel contact with one of 
its counterparts and only a capacitive coupling with the other. From the point of view 
of Kondo effect such configuration can be treated as an extension of a T-shape quantum 
dot |H3 EHl IZni mi EH] ■ Triangular trimers of Cr ions on a gold surface were also studied 
PH| inSl inn IHH] • The electronic and magnetic structure of these trimers is described in 
terms of a three-site Kondo effect. The orbital symmetry of triangle is discrete. It results 
in additional degeneracies of the spectrum of trimer, which may be the source of non- 
Fermi-liquid fixed point PH} lUUj. 

In this Chapter we concentrate on the point symmetry of ring-like TTQD and its 
interplay with the spin rotation symmetry in a context of Kondo tunneling through this 
artificial molecule. Indeed, the generic feature of Kondo effect is the involvement of 
internal degrees of freedom of localized "scatterer" in the interaction with continuum 
of electron-hole pair excitations in the Fermi sea of conduction electrons. These are 
spin degrees of freedom in conventional Kondo effect, although in some cases the role of 
pseudospin may be played by configuration quantum numbers, like in two-level systems 
and related objects 183]. TTQD may be considered as a specific Kondo object, where 
both spin and configuration (orbital) excitations are involved in cotunneling on an equal 
footing [TTinirTTTT] . 

To demonstrate this interplay, we consider a fully symmetric TTQD consisting of three 
identical puddles with the same individual properties (energy levels and Coulomb blockade 
parameters) and inter-dot coupling (tunnel amplitudes and electrostatic interaction). Like 
in the above mentioned triangular ratchet IM] . we assume that the TTQD in the 
ground state is occupied by one electron and Coulomb blockade is strong enough to 
completely suppress double occupancy of any valley j = 1, 2, 3. This means that the only 
mechanism of electron transfer through TTQD is cotunneling, where one electron leaves 
the valley j into the metallic leads, whereas another electron tunnels from the reservoir 
to the same valley j or to another valley /. In the former case only the spin reversal is 
possible, whereas in the latter case not only the spin is affected but also the TTQD is 
effectively "rotated" either clockwise or ant i- clockwise (see Fig. 15. Ij) . 

Discrete rotation in real space and continuous rotation in spin space may be encoded 
in terms of group theory. The group C^^ characterizes the symmetry of a triangle, and 
the group SU{2) describes the spin symmetry. As a result, the total symmetry of TTQD 
is determined by the direct product of these two groups. One may use an equivalent 
language of permutation group P3 for description of the configurations of TTQD with an 



67 




Figure 5.1: Clockwise (c) and anti-clockwise (a) "rotation" of TTQD due to cotunneling 
through the channels 2 and 3, respectively. 

electron occupying one of three possible positions in its wells. The discrete group P3 is 
characterized by three representations A, B, E or three Young tableaux [3], [1'^], [21]. Here 
A and B are one- dimensional representations with symmetric and antisymmetric basis 
functions, respectively, and E is the two-dimensional representation with basis functions, 
which are symmetric or antisymmetric with respect to mirror reflections. In the configu- 
ration shown in Fig. 15.2b . the perfect triangular symmetry characterizes both the triple 
dot and three leads, whereas in the configuration illustrated by Fig. 15.2b . the permutation 
symmetry between the sites (2,3) and (1,3) is broken, while the system remains invariant 
relative to (1,2) permutation. 




(a) (b) 

Figure 5.2: Triangular triple quantum dot (TTQD) in three-terminal (a) and two-terminal 
(b) configurations. 

Since the TTQD has a ring geometry, it is sensitive to an external magnetic field 
directed perpendicular to its plane, because the electron acquires additional gauge phase 
as it tunnels between the wells. The phase is determined by the magnetic flux through the 
ring. In a two-terminal geometry (Fig. 15.2b ). the Aharonov-Bohm (AB) effect is possible 
due to the interference between the channels (sl3d) and (s23d). Since the magnetic 
flux in this case is separated between two rings, the AB oscillations should have more 
complicated character than in a text-book AB interferometer. The U{1) symmetry of 
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an electron in an external magnetic field also influences the Kondo resonance because it 
changes the total symmetry of TTQD (breaks the chiral symmetry). In this Chapter we 
show that the interplay between continuous SU (2) symmetry, discrete C^^ symmetry and 
gauge U{1) symmetry can be described in terms of dynamical symmetry group approach 
ISa ESI EZl ESI im discussed in Chapters 3 and 4. 

It should be stressed that the interplay between the Kondo tunneling and AB inter- 
ference differs from the effects considered recently in the context of mesoscopic quantum 
interferometer containing a Kondo impurity on one of its arms jlU2[ I1U3| llU4t llU5t llU(jj . 
Unlike the Fano-like effects which take place in the latter case, the coherent TTQD as a 
whole plays part in the physics of an AB interferometer and a resonance Kondo-scatterer. 

5.2 Hamiltonian 

A symmetric TTQD in a contact with source and drain leads (Fig. 15. 2j) is described by 
the Anderson Hamiltonian for lead electrons Ckja and dot electrons dja, 

H = H^ + Hi,ad + Ht. (5.1) 

The flrst term, H^, is the Hamiltonian of the isolated TTQD, 

3 

3 = 1 <T j (jl) <T (jl) cr 

where o =t, i is the spin index, (jl) = (12), (23), (31). Here Q and Q' are intra-dot 
and inter-dot Coulomb blockade parameters {Q ^ Q'), and W is the inter-dot tunneling 
parameter. The second term, Hi^ad, describes the electrons in the leads labelled by the 
same indices j = 1,2,3 as the dots in the case of Fig. 15.2k . In the case of Fig. 15.2b 
j = s,d for source and drain electrodes respectively, 

Hlead = '^'^^kjclj^Ckja- (5.3) 

fccr j 

The last term, Ht, is the tunneling Hamiltonian between the dot and the leads. It has 
the form 

Ht = {^A^A^ + H.c) , (5.4) 

fccr j = l,2,3 

in the 3-terminal geometry (Fig. 15.2b ) and 

Ht = {VsAsA<^ + vAdJ^^r + H.c) (5.5) 

fccr i = l,2 
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in the 2-terminal geometry (Fig. 15.2b ). We assume that in the latter case the mean-free 
path for electrons near the "tip" of the source electrode exceeds the size of this tip. The 
2-terminal device has the symmetry of isosceles triangle with one mirror reflection axis 
(1^2). 

First we consider TTQD with three leads (Fig. 15. 2h .) whose ground state corresponds 
to a single electron occupation A/" = 1, and assume that all three channels are equivalent 
with V <^ W so that the tunnel contact preserves the rotational symmetry of TTQD, 
which is thereby imposed on the itinerant electrons in the leads. It is useful to re-write the 
Hamiltonian in the special basis which respects the C^^ symmetry, employing an approach 
widely used in the theory of Kondo effect in bulk metals |107t llU8j . The Hamiltonian 
Hd + Hiead is diagonal in the basis 

4,. = ^ {dl + dl + dl) , 4^^^ = i= (dl + e^'^^dl + e^^^dl) ; (5.6) 

Here = 27r/3, while A and E form bases for two irreducible representations of the 
group Cs^. Only a symmetric representation A of the P3 group arises in this case. The 
antisymmetric state B cannot be constructed due to the well known frustration property 
of triangular cells. The spin states with A/" = 1 are spin doublets (D), so the Hamiltonian 
of the isolated TTQD in this charge sector has six eigenstates. They correspond to a 
spin doublet \A) with fully symmetric "orbital" wave function {A) and two degenerate 
doublets \E±). The corresponding single electron energies are, 

EDA = e + 2W, EDE = e- W. (5.8) 

The Anderson Hamiltonian ()5.1|1 rewritten in the variables ()5.6|1 and ()5.7j) may be 
expressed by means of Hubbard operators X^^' = |A)(A'|, with A = 0, F, A: 

H = j2 Exx'' + ^^^r,k + Yl [y^'x^'^^M + E ^'^4',..^'^ + h.c] . (5.9) 

A ka r,k r,kcr AV 

Here |0) stands for an empty TTQD, |F) = \DA), \DE) belong to the single electron 
charge sector, and |A) are the eigenvectors of two-electron states. The eigenstates E'a for 
A/" = 2 are 

EsA = e2 + 2W ETE = e2 + W, 

EsE = e2-W ETB = e2-2W. (5.10) 

Here €2 = 2e + Q', the indices S,T denote spin singlet and spin triplet configurations of 
two electrons in TTQD, and the inequality W Q is used explicitly. The irreducible 



70 



representation B contains two-electron eigenfunction, which is odd with respect to per- 
mutations j ^ I. The tunnel matrix elements are redefined accordingly, 

A peculiar feature of ring configuration is an explicit dependence of the order of levels 
within the multiplets ()5.8j) and ()5.1()|1 on the sign of the tunnel integral W . For M = 1 
the doublet \DA) is a ground state, provided W < In the case of > the lowest 
levels are the doublets \DE±). The orbital degeneracy of the states E± is a manifestation 
of rotation/permutation degrees of freedom of the TTQD. In the next section we will 
show that these discrete rotations are explicitly involved in Kondo tunneling. 



5.3 Magnetically Tunable Spin and Orbital Kondo 
Effect 

To describe the infiuence of an external magnetic field B (perpendicular to the TTQD 
plane) on the dot spectrum, one may treat the TTQD as a three-site cyclic chain with 
nearest-neighbor hopping integrals — connecting these sites. The spectrum of this 
"chain" is 

EDr{p) = e-2Wcosp, p = 0, 2n/3, An/S. (5.11) 

This equation is the same as ()5.8|) with negative W, where p = 0, 27r/3, An/ 3 correspond 
respectively to F = A,E+,E^. A perpendicular magnetic field modifies this spectrum. 
If the Zeeman splitting is weak, the only effect of the magnetic field is refiected by an 
additional phase acquired by the tunneling integral W. This phase is determined by the 
magnetic flux $ through the triangle, so that the spectrum of the TTQD becomes, 

Eorip,^) =e-2Wcos(^p-^y (5.12) 

Figure E31 illustrates the evolution of Ejjrip, $) induced by B. Variation of B between 
zero and Bq (the value of Bq corresponding to the quantum of magnetic flux $o through 
the triangle) results in multiple crossing of the levels -E^r- The periodicity in $ of the 
energy spectrum is 2tt since each crossing point (Fig. 15. Hj) may be considered as an orbital 
doublet E± by regauging phases ip in Eq. ()5.6|) . 

The accidental degeneracy of spin states induced by the magnetic phase $ introduces 
new features into the Kondo effect. In the conventional Kondo problem, the effective 
low-energy exchange Hamiltonian has the form JS ■ s, where S and s are the spin op- 
erators for the dot and lead electrons, respectively Here, however, the low-energy 
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states of TTQD form a multiplet characterized by both spin and orbital quantum num- 
bers. The effective exchange interaction reflects the dynamical symmetry of the TTQD 
[S3 ESI EH 113 mi • As was mentioned in Sec. 12.41 the corresponding dynamical symmetry 
group is identified not only by the operators which commute with the Hamiltonian but 
also by operators inducing transitions between different states of its multiplets. Hence, 
it is determined by the set of dot energy levels which reside within a given energy in- 
terval (its width is related to the Kondo temperature Tk). Since the position of these 
levels is controlled in this case by the magnetic field, we arrive at a remarkable scenario: 
Variation of a magnetic field determines the dynamical symmetry of the tunneling de- 
vice. Generically, the dynamical symmetry group which describes all possible transitions 
within the set {DA, DE-i-} is SU{6). However, this symmetry is exposed at too high 
energy scale ~ W, while only the low-energy excitations at energy scale <^ W are 
involved in Kondo tunneling. It is seen from Fig. 15. HI that the orbital degrees of freedom 
are mostly quenched, but the ground state becomes doubly degenerate both in spin and 
orbital channels around $ = {2n + 1)tt, {n = 0, ±1, . . .). 




2.5 5 7.5 10 12.5 15 17.5 

Figure 5.3: Upper panel: Evolution of the energy levels Ea (solid line) and E± (dashed 
and dash-dotted line, respectively.) Lower panel: corresponding evolution of Kondo tem- 
perature. 

Let us compare two cases: $ = 0, (the ground state is a spin doublet DA), and 
$ = TT, (the ground state is both orbital and spin doublet DE). It is useful at this 
point to generalize the notion of localized spin operator S"* = \a)Ti{a'\ (employing Pauli 



72 



matrices fj {i = x, y, z)) to Sfp, = \Ta)fi{a'T'\ (F, V = A, E±), in terms of the eigenvectors 
()5.6|) . Similar generalization applies for the spin operators of the lead electrons: Spp/ = 
J2kk' ka^i^^',k'a' with Cr^ka defined in ()5.7|) . First, when $ = 0, the rotation degrees 
of freedom are quenched at the low-energy scale. The only vector, which is involved in 
Kondo co-tunneling through TTQD is the spin Saa- Applying SW procedure jSH], the 
effective exchange Hamiltonian reads, 

Hsw = Je (Saa ■ SE_f_E+ + Saa ■ se_e.) + -'^aSaa ■ saa- (5.13) 

The exchange vertices Jr are 

2V^ f 1 1 \ , , 



3 yep - e e + Q-ep e + Q' - ep , 
Note that > as in the conventional SW transformation of the Anderson Hamiltonian. 
On the other hand, Jp < due to the inequality Q ^ Q'. Thus, two out of three available 
exchange channels in the Hamiltonian ()5.13|) are irrelevant. As a result, the conventional 
Kondo regime emerges with the doublet DA channel and a Kondo temperature, 

T^^^ = Dexp|-^|. (5.15) 

where ja = PqJa- 

Second, when $ = vr, the doublet DA is quenched at low energy, and the Kondo effect is 
governed by tunneling through the TTQD in the state \DE) whose symmetry is SU{A). 
This scenario of orbital degeneracy is different from that of occupation degeneracy studied 
in double quantum dot systems |l()9j . The 15 generators of SU (4) include four spin vector 
operators S^;^^;;^ with a, 6 = ± and one pseudospin vector T defined as 

T+ = Y,\E+,cr){a,E^\, = ^ (|i5;+, E+| - i?„|) . (5.16) 

Its counterpart for the lead electrons is 

a cr 

Due to SU{4:) symmetry of the ground state, the SW Hamiltonian acquires a rather 
rich structure, 

Hsw = Jl{^E+E+ ■ ^E+E+ + Se-E- ■ ^E-E-) 

+ <^2(Se+E+ ■ ^E-E- + Sp_p_ ■ Sp_^_E^) + J3{Se+E+ + ^E-E-) ' ^AA 

+ J4{Se+E^ ■ Se^E+ + Se^E+ ■ ^E+E^) (5.18) 

+ J5{^E+E- ■ {SaE. + Se+a) + Se.E+ " {SaE+ + Se_a)) + JqT ■ T, 
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where the couphng constants are 

1/2 

Jl = Ji = JA, J2 = J-i = J^ = Je, J& = \ • (5.19) 

Thus, spin and orbital degrees of freedom of TTQD interlace in the exchange terms. The 
indirect exchange coupling constants arise due to co-tunneling processes with virtual exci- 
tations of states with zero and two electrons. These constants include both diagonal (jj) 
and non-diagonal (j7) terms describing reflection and transmission co-tunneling ampli- 
tudes. (Our representation of spin operators and therefore the form of spin Hamiltonians 
flOT?|) . (fOT^ differs from that used in [HSlinn] and in pimmi] .! 

The interplay between spin and pseudospin channels naturally affects the scaling equa- 
tions obtained within the framework of Anderson's "poor man scaling" procedure |5f)j . 
The system of scaling equations has the following form: 



dt 



,■2 a2 



dt 



r 2 2~l r / \ / \i 

= - [h + h\ ^ = ~ ^^(ji + J2 + Je) + Jeiji - J2)\ , 

Here ji = poJi, and the scaling variable is t = In(po-D). Analysis of solutions of the 
scahng equations ()5.20|) with initial values of coupling parameters listed in Eq. ()5.19|) . 
shows that the symmetry-breaking vertices js and js are irrelevant, but the vertex j2, 
which is negative at the boundary D = D evolves into positive domain and eventually 
enters the Kondo temperature. The latter is given by the following equation 

) = Dexp I ] . (5.21) 

We see from 1)5.211) that both spin and pseudospin exchange constants contribute on an 
equal footing. Unlike the isotropic Kondo Hamiltonian for A/" = 3 discussed in Refs. 
IHHl innj, the non- Fermi-liquid regime is not realized for A/" = 1 with Hsw ()5.18|) . The 
reason of this difference is that starting with the Anderson Hamiltonian ()5.1|) with finite 
Q, Q', one inevitably arrives at the anisotropic SW exchange Hamiltonian for any A/". As 
a result, two out of three orbital channels become irrelevant. However is enhanced due 
to inclusion of orbital degrees of freedom. Moreover, this enhancement is magnetically 
tunable. It follows from ()5.12|) that the crossover SU (2) SU{A) SU (2) occurs three 
times within the interval < $ < Gvr and each level crossing results in enhancement of 
Tx from ()5.15|) to 1)5.211) and back (lower panel of Fig. 15.3)1 . These field induced effects 
may be observed by measuring the two-terminal conductance Gji through TTQD (the 
third contact is assumed to be passive). Calculation by means of Keldysh technique (at 
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T > Tk) similar to that of Ref. [58] show sharp maxima in G as a function of magnetic 
field, following the maxima of Tk (Fig. I5.4|l . 

G/Go 
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Figure 5.4: Evolution of conductance {Go = ne'^/h). 

So far we have studied the influence of the magnetic field on the ground-state symmetry 
of the TTQD. In a two-lead geometry (Fig. 15.2b ) the field B affects the lead-dot hopping 
phases thereby inducing the AB effect. (The necessary condition for AB effect is that the 
electron coherence length in the source should exceed the size of the electrode "tip".) The 
symmetry of the device is thereby reduced since it looses two out of three mirror reflection 
axes. The orbital doublet E splits into two states, but still, the ground state is \DA). In 
a generic situation, the total magnetic flux is the sum of two components $ = $i -|- $2- 
In the chosen gauge, the hopping integrals in Eqs. ()5.2|) . ()5.5|) are modified as, 

W ^ Vrexp(2$i/3), ^1,2 ^ Kexp[±^($i/6 + $2/2)], 

and the exchange Hamiltonian now reads, 

H = JsS ■ Ss + JdS ■ Sd + JsdS ■ (Ssd + Sds). 



(5.22) 



Here 



J. 



Jd 



/l + cos + $2) 



(-F 



+ 



t + Q ~ ep 



Jsd 

for < $1 < TT, and 



^VsVd 



1 



e + Q - ep 
1 



(5.23) 



+ 



e + Q'-ep 



cos 



$1 

6 



Js 
Jd 
Jsd 



If 1^ 1 - sin $2 
3 



4Vf 



3 \ep — e 
2VsVd ( 1 



e + - cf e + Q - eF 
1 1 



< 0, 



+ 



t^Q - ep 
1 



(5.24) 



+ 



t^Q' -tp 



sin $2, 
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for $1 = 71. Applying poor man scaling RG procedure on the Hamiltonian ()5.22j) yields 
the scahng equations, 



djs 



dlnD 

djd 
dlnD 

djsd 



dlnD 

Eqs. fj5.25|) give the Kondo temperature, 



- [fs+Jd] , 

~ [jd + Jd] 5 

-jsd \js + jd] ■ 



(5.25) 



Tk = Do exp 



The conductance at T > reads 

G_ 



Js + 3d + ^{]s-]dY+^]l 



(5.26) 



ild 



(5.27) 



A{3s+3dYln\T/TK)' 

It was verified that the resulting conductance (?($!, $2) ()5.27p obeys the Byers-Yang 
theorem (periodicity in each phase) and the Onsager condition G'($i, $2) = G'(— $1, — $2)- 
We choose to display the conductance along two lines <l>i(<l>), $2($) in parameter space of 
phases, namely, G($i = $, $2 = 0) and G{^i = <l>/2, $2 = $/2) (figure l531 left and right 
panels, respectively). The Byers-Yang relation implies respective periods of 2tt and An in 
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Figure 5.5: Conductance as a function of magnetic field for $2 
$1 = $2 = $/2 (right panel). 



(left panel) and 



$. Experimentally, the magnetic fiux is of course applied on the whole sample as in figure 
lb, and the ratio $i/$2 is determined by the specific geometry. Strictly speaking then, 
the conductance is not periodic in the magnetic field unless $1 and $2 are commensurate. 

The shapes of the conductance curves presented here are distinct from those pertaining 
to a mesoscopic AB interferometer with a single correlated QD and a conducting channel 
[TIM ITM ITT^ (termed as Fano-Kondo effect [Ml). For example, G($) in Fig. 3 of 
Ref. |l()3j (calculated in the strong coupling regime) has a broad peak at $ = 7r/2 with 
G($ = vr/2) = 1. On the other hand, G'($) displayed in Fig. 15.31 (pertinent to Fig. 
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15.2b and obtained in the weak coupling regime), is virtually flux independent except near 
the points $ = {2n + l)7r (n integer) at which the SU{4:) symmetry is realized and G is 
sharply peaked. The phase dependence is governed here by interference effects on the level 
spectrum of the TTQD. The three dots share an electron in a coherent state controlled by 
the strong correlations, and this coherent TTQD as a whole is a vital component of the AB 
interferometer. As a result, in the setup of Fig. 15.2b . the conductance vanishes identically 
on the curve Jgdi^i, ^2) =0. The AB oscillations result of interference between 

the clockwise and anticlockwise "effective rotations" of TTQD in the tunneling through 
the {13} and {23} arms of the loop (Fig. 15.2b ). provided the dephasing in the leads 
does not destroy the coherence of tunneling through the two source channels. (Interplay 
between two chiral states results in SU{4) Kondo effect in carbon nanotubes in axial 
magnetic field jll3t I114j .) On the other hand, in the calculations performed on Fano- 
Kondo interferometers, G{^) remains finite |l()3j . Incidentally, there should be a Fano 
effect due to the renormalization of electron spectrum in the leads induced by the lead-dot 
tunneling similar to that in chemisorbed atoms |115j but that is beyond the scope of this 
Thesis. 

5.4 Conclusions 

To conclude, we have shown that the Kondo and AB effects in TTQD expose peculiar 
symmetries in the physics of strongly correlated electrons. Dynamical symmetry is a 
universal tool, which allows to derive the effective spin Hamiltonians describing low-energy 
cotunneling with spin reversal through TTQD. This artificial molecule possesses both the 
continuous rotation symmetry in spin space and discrete rotation symmetry in real space. 
This discrete symmetry is imposed on the whole system both in the case of molecular 
trimer Crs on the metallic surface and in the three-terminal device - equilateral TTQD 
in its center (Figs. 15.11 15.2b ). The tunnel problem is mapped on the Coqblin-Schrieffer 
(CS) model of magnetic impurity with "orbital" degrees of freedom |l()7[ll()8] . Like in the 
latter case, the symmetry of the Kondo center is SU{2n). The parameter n characterizes 
additional orbital degeneracy in the conventional CS model. It describes the pseudospin 
operator of finite clockwise/anti-clockwise rotations of TTQD in the process of electron 
cotunneling. 

Abstract concepts like dynamical symmetries are shown to be directly related to exper- 
imentally tunable parameters. The conductance is sensitive to an external magnetic field 
both in three- and two-terminal geometry even in the weak coupling regime at T > Tk- 
In the former case, the conductance can be enhanced due to change of the dynamical 
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symmetry caused by field-induced level crossing (Fig. 15. 3|) . In the latter case, the conduc- 
tance can be completely suppressed due to destructive AB interference in source-drain 
cotunneling amplitude (Fig. 15. 5|) . These results promise an interesting physics at the 
strong coupling regime as well as in cases of doubly and triply occupied TTQD. It would 
also be interesting to generalize the present theory for quadratic QD which possesses 
rich energy spectrum with multiple accidental degeneracies. 
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Appendix A 

Diagonalization of the Trimer 
Hamiltonian 



Here we describe the diagonalization procedure of the Hamiltonian of the isolated TQDs 
occupied by four and three electrons (Fig. I4.2|) . The dot Hamiltonian has the form, 

= J2 Yl ^^^i-daa + J2 '^a^aT^ai + J] i^Ajaa + H.C.). (A.l) 

a=l,r,c (J a a=l,r 

(a) Four electron occupation: 

The Hamiltonian (jA.ljl can be diagonalized by using the basis of four-electron wave func- 
tions 

\ta, 1) = di^di^dyiio), 1) = d^didyiio), 

\^a) = ^{dl,dl-dldl) 4,4J0), 

lex) = 444^4J0), (A.2) 

where a = l,r] I = r, r = I. The Coulomb interaction quenches the states with two 
electrons in the central dot and we do not take them into account. The states ()A.2j) form 
a basis of two triplet and three singlet states. In this basis, the Hamiltonian ()A.1|) is 
decomposed into triplet and singlet matrices, 

^-{::} 

and 

/ El V2Wi \ 

Hs= Br \^Wr , (A.4) 
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where ei = ec + ei + le^ + Qr, £r = + e-r + 2e; + Qi, and e^e^ = 26^ + 2e,. + + Qr- 
We are interested in the hmit /?a -C 1 = Wa/{ea + Qa — Cc))- So the secular matrix 
may be diagonahzed in lowest order of perturbation theory in The eigenfunctions 
corresponding to the energy levels ()3.1|) are, 

\Sa) = V^r^2;52 \Sa) - V2f3a\ex), 

\Ta) = \ta) , (A.5) 

\Ex) = ^1-2/32 -2/32 I ex) + V2(3i\si) + V2(3r\sr). 

In the completely symmetric case, ei = Er = e, Qi = Qj. = Q, Wi = Wr = W, the 
eigenfunctions corresponding to the energies ()3.2|) are 

\Sl) + \Sr) 



\S^) = v/r^^^^^-2/3|ex), 

|^_) = M^fll, (A.6) 



\T±) 



V2 

\tl) ± \tr) 



V2 ' 

\Ex) = ^/l-4f3\ex) + V2(3{\si) + 

where /3 = iy/(e + Q-ec). 
(b) Three electron occupation: 

In this case the Hamiltonian (jA.lll can be diagonalized by using the basis of three-electron 
wave functions 

„ , ([rf+4 - rf+4]ci+ + - rf,v+]Oio) 

Ka) = , 

1^2,^) = - d+ |0), 

g,±^^=ci+±rf,V+|0), 



2, 



2/ 

a) = |0), a) = ci+ rf+ c/+ |0), 

a) = rf+ d+ d+ 10), |6„ a) = 4rf+c/+ |0), (A.7) 

where a =!> i- The three-electron states |A) of the TQD are classified as a ground state 
doublet l-Bi), low-lying doublet I-B2) and quartet \Q) excitations, and four charge-transfer 
excitonic doublets Bac and Ba (a = Z,r). In the framework of second order perturbation 
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theory with respect to jSa ()2.58|) . the energy levels E/^ are 

Eb, =ec + ei + er-^ [Wi(3i + W^f^r] , 

Eb^ = ec + ei + er-^ [Wipi + WrPr] , 
Eq = ec + ei + er, 

Eb^, = Cc + 2ea + Qa- W-a^a, 

Eb^ =ea + 2e-a + + WaPa + ^WajSa. (A.8) 

The eigenf unctions corresponding to the energy levels ()A.8|1 are the following combi- 
nations, 



|5i,cr) = 


lilbi, 


, \/6 ^ , \/6 ^ 

- ^A|Or,(T) + —/3r\bi,a^ 


|i?2,a) = 


72 1^2, 


(^) - ^Pi\br,cr) - ^pr\bi,(T[ 


\Q,s,) = 






\Bac, cr) = 






\Br,a) = 




2Pf-P?\br,a) + Pr\bic,cr) + - 


\Bi,a) = 




2(3^-(3f\bi,a)+(3i\br,c,a)- 



(A.9) 



V2, 



where 71 and 72 are determined by Eg. ()4.2H) . 
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Appendix B 

Rotations in the Source-Drain and 
Left-Right Spaces 



In the generic case, the transformation which diagonahzes the tunnehng Hamiltonian ()3.4p 
has the form 



QofcfT 



( 



\ / Q,fe, \ 



V 



Cldka 
y Crdka j 



(B.l) 



Ul Vi 

-vi ui 

Ur Vr 

—Vr Ur j 

with Ua = Vas/Va, fa = Kd/K; = + = l,r).liaa symmetric case = = 
V, this transformation simphfies to 

(^aeka 2 ^ i^aska ~l~ (^adku) i (^aokcr 2 ( Cdgj^^j -\- Cg^dka^ i (-B'^) 



and only the even (e) combination survives in the tunnehng Hamiltonian 

Ht^rr. = VY,icLkJaa + H.C.). 



(B.3) 



aka 



So the odd combination (o) may be omitted. 

If, moreover, the whole system "TQD plus leads" possesses / — r symmetry, ei = Er, 
the second rotation in / — r-space 



(B.4) 





^gka 






/ 


1 


1 








\ 










C-uka 




1 




-1 


1 












C-rekcr 






dga 




"7! 










1 


1 






dlcr 




v 


dufj 


I 




V 








-1 


1 


/ 


v 


dra 


1 



transforms Hiead + Htun into 

Hlead + Htun = ^ |^efe^n^fccr + ^(c|,fco-'^»?o- + H.C.] 

•qka 

with Cfcg = Cfc — tir, tku = Cfe + tlr- 



(B.5) 
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Appendix C 

Effective Spin Hamiltonian 



The spin Hamiltonian of the TQD with A/" = 4 occupation in series geometry fFig J4.2|) is 
derived below. The system is described by the Hamiltonian (jOJ. The Schrieffer- Wolff 
transformation [221 for the configuration of four electron states of the TQD projects out 
three electron states |A) and maps the Hamiltonian ()4.H) onto an effective spin Hamilto- 
nian H acting in a subspace of four-electron configurations (A| . . . |A'), 

H = e'^He-'^ = + ^ [S..\S, H]]..], (C.l) 

m 

where 

T/AA 

5 = V V X^\,,, + H.c. (C.2) 

AA (k).,a - ^^'^ 

Here (k) stands for the electron or hole states whose energies are secluded within a layer 
±D around the Fermi level. E\x = E\{D) — E\{D) and the notation a = /,r is used. 
The effective Hamiltonian with three-electron states |A) frozen out can be obtained by 
retaining the terms to order 0(|V"p) on the right-hand side of Eq. ljdlll . It has the 
following form, 

H = Y, EaX^"" + E ''^'^^k.cak. (C.3) 

A {k)a,a 

- X] X] X] Jkk^a^^^ clku^^ak'a' - X] X] i-Jkk^lrX^^ C^ka'^lk'a' + H.C.), 



where 



AA'A kk'aa' a=Lr AA'A kk'aa' 



jAA' _ f 1 ^ 1 

"^kk'a — n ' 



EaX — ^ka -E'a'A — ^k'a 

,-.(ffl(^,^). ,c.4, 

\^AX — ^kl -C/A'A — efc'r / 

The constraint YIa-^^^ = 1 is valid. Unlike the conventional case of doublet spin 1/2 we 
have here an octet A = {A;, A,.} = {Si, Ti, Sr, T^}, and the SW transformation intermixes 
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all these states. The effective spin Hamiltonian ()C.3|) to order 0(|l^p) acquires the form 
of Eq.dOJ. 
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Appendix D 
SO (7) Symmetry 



The operators S;, S^, R^, Ri, R2, R3 and Ai {i = 1,2,3) (see Eas. dTTHll . innnil . dCTI ^ 
obey the commutation relations of the 07 Lie algebra, 

[Rlji Slk] = iejkmRlm, [Rlj, Srk] = [Rsj, Sik] = 0, 

[-Rsj; -Rsfc] '^(ijkmSrmi \R'iji ^Tk\ ^(^jkm-Rsm^ 

i 1 

[Rlj, Rlk] = iejkmSrmi'^ — Sjz)i^ — 5kz) + -ejkmSlm{Sjz + Skz) — -{SljSkz — SlkSjz), 

i 1 

[i?2j; R2k] = i^jkmSlmiX ~ ^jz)i)- ~ ^kz) + '^^jkmSrmi^jz + ^kz) ~ '^{S^jSkz ~ S^-k^jz), 

i 1 
[Rlj, R2k] = 2^jkm{SrmSjz + Sim^kz) + [Slj^kz ~ Srk^jz + (Siz — Srz)Sjz5kz] ; 

[Rsj, Rlk] = i^jkmRlmi^ — Sjz —) — 5jz)Rlj, 

[Rsj, R2k] = i^jkmRlmiSjz + + ^ Sjz)Rlj, 

[Rlj, Rlk] = i^jkmRsmiSkz + ^(1 — ^kz)Rzk, 

[R2j, Rlk] = i^jkmRzm{^ — ^kz ^) + 'l^i^ ~ ^kz)Rzk, 



[Ai, Sij] = iA26jz + -^{RixSjx - Rly^, 



iy/2 ~ 

[A2, Sij] = —iAi5jz ^{Riy^jx + RixSjy), 

[Al, Srj] = iA25jz —{R2x^jx ~ R2ySjy), 

[A2, Srj\ = —iAiSjz H —{R2ySjx + R2xSjy), 

[A3, Sij] = -iR2j{l - 6jz), [A3, Srj] = iRij{l - Sjz), 
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[Ai, Ri^ = —{Rsx^jx — Rsy^jy), [-^2, Rlj] = —^{Rsy^jx + R^x^jy)-, 

[^1, -Rsj] = —^{Rlx^jx — Rly^jy), [^2, Rsj] — ^{Rly^jx + Rlx^jy), 

[A3, Rlj] = —iRszSjz, [■^3,R3j] = iRlzSjz, 

[Ai, Rlj] — ——^{SixSjx — SiySjy), [A2, R-ij] — —^{SiySjx + SlxSjy), 

[A3, Rlj] = ~i{SrxSjx + Sry5jy)i [^li R2j] — ^ {^rx^jx ~ Sry3jy)i 

[A2, R2]\ = ^{^rySjx + SrxSjy), [A^, i?2j] = i{SlxSjx + Sly5jy), 

[Ai, A2] = -t{Si, + Sr.), [Ai, A3] = [A^, A3] = 0, 

[Saj, R,k] = r;ZR,m + [^3,-, Rlk] = P^,mR,m + a],A^. (D.l) 

Here j,k,m are Cartesian indices, a = I, r; iJ,,^ = 1,2; n = 1,2,3; r"^^, a"^", a"^ and 
Pjkm are the structural constants, rj^;; = rj^, q;^^ = -a^^j^" (1 = 2, 2 = 1). Their 
non-zero components are: 
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Hyyz 
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" 2' 
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rxzy 


rzyx 




0^ -- 
ryzx 
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rzxy 





~ ^ZZ ~ ^V^, (ij-J^ — Cty^ — Otyy " Z . 
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The following relations hold, 



S„ • = S„ • R3 = 0, = = 0, 

S2 ^ 2X'^«^», Ri • R| + R2 • = 2 ^ 

a=/,r 

r2 = X''"'' + 3X^'^' , R| = X^'"''^ + 3X^'^' , 

+ + = + JS^'''■'^^ (D.2) 

Therefore, the vector operators S^, S^, R;, Ri and scalar operators {i — 1,2,3) generate 
the algebra 07 in a representation specified by the Casimir operator 

Sf + + Rf + ^ R, • Rl + R^ + J] = 6. (D.3) 
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Appendix E 

Young Tableaux Corresponding to 
Various Symmetries 

A TQD with "passive" central dot and "active" side dots reminds an artificial atom with 
inner core and external valence shell. The many-electron wave functions in this nano- 
object may be symmetrized in various ways, so that each spin state of M electrons in 
the TQD is characterized by its own symmetrization scheme. One may illustrate these 
schemes by means of the conventional graphic presentation of the permutation symmetry 
of multi-electron system employing Young tableau |117j . For instance, triplet state of two 
electrons which is symmetric with respect to the electron permutation is labelled by a 
row of two squares, whereas the singlet one which is antisymmetric with respect to the 
permutation is labelled by a column of two squares. Having this in mind we can represent 
the singlet and triplet four electron states of the TQD ()A.5|) by the four tableaux shown 
in Fig. IE. II The tableaux Si {Sr) and Ti {%) correspond to the singlet and triplet states 
in which the right (left) dot contains two electrons (grey column in Fig. IE.1|) whereas 
electrons in the left (right) and central dots form singlet and triplet, respectively. 



S, T, 



Figure E.l: Young tableaux corresponding to the singlet [So) and triplet (T^) four electron 
states of the TQD. The grey column denote two electrons in the same dot (right or left). 

The Young tableaux corresponding to various SO{n) symmetries discussed in Chapters 
3 and 4 can be obtained by combining the appropriate tableaux (Fig. IE.2|1 . The highest 
possible symmetry P x 50(4) x SO {A) is represented by four tableaux T;, T^, Si and Sr 
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Px50(4)x50(4) 

PxSO{3)xSO{3) 
SO{l) 




^0(5) ^ ^ 



50(4) 




Figure E.2: Young tableaux corresponding to SO{n) symmetries. 

since all singlet and triplet states are degenerate in this case. The symmetry P x SO{?>) x 
5*0(3) occurs when two triplets Ti and T,. are close in energy and these are represented 
by the couple of Young tableaux in the second line. Following this procedure, the 5*0(7) 
symmetry can be described in terms of two triplets Ti, diagrams and one singlet Si 
diagram. Moreover, 5*0(5) symmetry is represented by two singlet 5*;, 5*^ diagrams and 
one triplet Ti diagram and, finally, one triplet and one singlet tableaux correspond to the 
5*0(4) symmetry. 
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